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A Associação dos Estudantes 
do Instituto Superior Técnico, a AEIST, 
foi fundada em 11 de Dezembro de 1911 


No tempo, 75 anos que são passados, 

a fundação do Instituto e a fundação da Associação coincidem. 
Praticamente são factos que emergem na mesma data e na mesma atitute. 
Prenúncio de que o Instituto e a sua Associação de Estudantes 

ficariam ligados. 


De facto assim tem sido. 

A Associação nasceu do Técnico, 

mas o Técnico, o Instituto, deve parte importante da sua configuração e 
prestígio ao que a sua Associação de Estudantes tem sabido e 

desejado ser. 


Na sessão de 11 de Dezembro comemorou-se a fundação da Associação dos Estudantes do I.S.T. 
Este número da TECNICA é especialmente dedicado a essa efeméride. 
Quanto se disse nessa cerimónia aqui fica para sinal e reflexão. 


O belo retrato de Alfredo Bensaúde, do número anterior, parece que nos chega de um passado 
distante. 

Medem-se os 75 anos do I.ST. e da A.E.I.ST. no ar antigo desse retrato. 

No entanto, as notas Histórico-Pedagógicas são de modernidade que ainda não se cumpriu. 


Neste número, o diamante 

simbolo das Bodas de Diamante 

invade com a sua luz o futuro. 

Nesse futuro, como no passado, tudo continuará a ser trivial e simples: 


É estudante quem estudar. 

E professor quem ensinar. 

Que possa estudar quem estuda e que possa ensinar 
quem saiba e se ofereça. 


Que o Instituto Superior Técnico e a Associação dos Estudantes do I.ST., unidos e em 
preocupação comum, cuidem destes factos simples e que exactamente por o serem definirão o 
futuro. 


Manuel José de Abreu Faro 


Palavras proferidas na Sessão de encerramento das comemorações dos 75 anos 
da Associação de Estudantes do I.ST. por 


A. SIMÕES LOPES 
Vice-Reitor no exercício das funções de Reitor da UTL 


Direi que é para mim extremamente grato aproveitar a oportunidade para saudar os alunos do 
Instituto Superior Técnico através da sua Associação; para saudar essa mesma Associação que acaba 
de completar 75 anos de existência e para, em reflexão breve, formular votos de muito empenhamento 
e de grande sucesso a todos quantos vêm dedicando parte significativa das suas energias ao movimento 
associativo, com intenções de organizar os estudantes para os grandes projectos do futuro. 

Sendo professor de outra Instituição da Universidade, sempre tenho apreciado, afirmado e defen- 
dido do Instituto Superior Técnico a imagem de uma grande Escola da Universidade Técnica e do país, 
qualquer que seja o índice que se utilize na aferição da grandeza. Essa imagem é fruto de um longo 
e persistente trabalho de mais de 75 anos, em que se vem tornando hábito referenciar figuras-chaves 
do passado, desde Alfredo Bensaúde e Duarte Pacheco, professores prestigiosos que muito prestigia- 
ram a Escola, havendo certa inibição, em referir os presentes, com o receio natural de a algum esque- 
cer e com a certeza de que a muitos feririamos na modéstia que anda normalmente de par com as 
qualidades positivas do cientista, do universitário, do Homem. Mas são os professores os mais fre- 
quentemente referidos, por razões óbvias, sempre se secundarizando aparentemente os alunos, nesse 
papel de acumulação do prestígio da Instituição. 

Não me consta que tenha havido, até hoje, contestação a esta atitude; e louvável (e correcto) é 
que não exista, não só porque nos velhos Mestres (e nos bons Mestres do presente) os alunos sempre 
sabem distinguir a imagem séria e digna dos que trabalham para a causa institucional, mas também 
porque neles se perfilam modelos de estudante, que o têm de ser sempre, porque a qualidade, a serie- 
dade, a dignidade da docência exigirão em pleno, em permanência, a preparação, o trabalho, o studio. 
Não pode conceber-se, muito menos aceitar-se, o professor que não estuda; o professor é antes do mais, 
e sempre, um estudante, até para merecer a consideração e o apreço dos alumni que lhe são confiados 
e se lhe confiam. 

Mas as reflexões que quero aproveitar para fazer vão noutra direcção, e penso que se justificam 
sem reservas. Quero deixá-las à consideração de todos, neste dia que é de festa e de celebração, mas 
em que olhar para o futuro não é menos importante do que glorificar o passado. 

Estará perfeitamente criado um verdadeiro espírito de comunidade, nesta grande Escola de cerca 
de 6000 alunos? Estaremos todos a dar o nosso melhor esforço no sentido de instalar, em bases perfei- 
tamente sólidas, um amplo espírito de comunidade nesta nossa Universidade onde o número de alunos 
anda perto dos 15000? E será importante criar-se esse espírito comunitário? 

Tenho-me dado conta, mais vezes do que desejaria, de que uma certa compartimentação existe 
entre os alunos desta Escola, mais do que em qualquer outra da Universidade, o que de certo, tem 
a ver com o tamanho da Instituição. O exemplo mais frisante surge-me a nível dos finalistas, sobre- 
tudo, quando começam a organizar as suas viagens de fim de curso para que se vai tornando habitual 
dirigirem-se à Reitoria para apoio — tantas vezes, infelizmente, para receberem apoio ''moral” ape- 
nas; e a observação que faço, regularmente, é a de que cada Departamento elabora o seu programa 
com grande alheamento do programa dos colegas de outros Departamentos, por vezes, mesmo, com 
largo desconhecimento deles. 

Direi que não é esse alheamento causa do frouxo apoio que a Reitoria pode dar; isto é, o apoio 
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talvez não fosse significativamente diferente se alguma coordenação houvesse entre os programas ela- 
borados. Mas, a mim, causa-me apreensão sincera o facto de a compartimentação existir, porque isso 
quer dizer que o convívio é menos intenso e diversificado do que eu desejaria (e se desejaria). Donde, 
o não aproveitamento de todas as potencialidades formativas que na Escola existem, porque ao enge- 
nheiro mecânico de amanhã (de hoje) há-de em muitas circunstâncias assistir a necessidade do traba- 
lho conjunto com o engenheiro químico, ou de minas,... que sei eu... amanhã, hoje mesmo... e, no 
entanto, virá a ser depois da Escola que há-de vencer-se a barreira para a actividade interdisciplinar, 
designadamente para o conhecimento dos Colegas dos Grupos que hão-de constituir-se. 

E a nível da Universidade? Salvo casos que me recuso mesmo assim a não considerar esporádicos, 
até para minha defesa do conceito de Universidade e da sua operacionalidade, que conhecimento alar- 
gado terão os alunos do “Técnico” dos seus colegas, e do trabalho dos seus colegas, das outras Esco- 
las? Será, por exemplo, que todos estarão seguros de conhecer quantas e quais são as outras Escolas 
da sua Universidade? 

E aqui também vem instalar-se a angústia, porque de novo se deixa de aproveitar as potencialida- 
des formativas existentes na Universidade, pela compartição excessiva que se instalou. Quantos alunos 
finalistas de engenharia civil, interessados no planeamento urbano, têm convivido com colegas seus 
da Faculdade de Arquitectura, do Instituto Superior de Agronomia (Arquitectura paisagista), do Insti- 
tuto Superior de Economia (Economia Regional e Economia Urbano), do Instituto Superior de Ciên- 
cias Sociais e Políticas (Sociologia Urbana e Sociologia Rural), interessados na mesma temática? 

E, no entanto, a actividade de cada um é cada vez menos destacável, cada vez menos isolável da 
actividade dos outros... 

Entretanto, falamos da necessidade de alargamento do conhecimento além-fronteiras, e defendêmo-lo, 
designadamente no sentido de aumentar a mobilidade dos estudantes que lhes permita viver e estudar 
noutros países da Comunidade Económica Europeia, como base (a mais sólida) de formação da Co- 
munidade dos Cidadãos. Mas estará feito o essencial para nos conhecermos e organizarmos, em Portu- 
gal, como Comunidade? 

Deixo estas reflexões a todos vós, porque embora sabendo que não é de vós que se espera a contri- 
buição total para a solução do problema, a verdade é que ela em parte de vós que se espera a contri- 
buição total para a solução do problema, a verdade é que ela em parte de vós depende, se a sensibiliza- 
ção existir e se a motivação despertar para um trabalho mais plural, mais conjunto, mais interligado. 
Sei que a Associação de Estudantes se empenha nesta perspectiva, em que tem actividade válida; e 
a testá-lo está a eleição recente para a Direcção da Associação Académica de Lisboa, que muito gosto- 
samente justificou de mim os cumprimentos e os desejos de bom labor; mas será esta minha preocupa- 
ção, esta preocupação da Direcção da vossa Associação, preocupação de todos? 

Retomo, para terminar, as minhas considerações iniciais. 

Encerrou-se, há semanas, o Ciclo de Comemorações dos 75 anos do Instituto Superior Técnico: 
foi motivo de alegria e de manifestação convicta de esperança no futuro desta grande Instituição. Ter- 
mina hoje o Ciclo de Comemorações dos 75 anos da Associação de Estudantes do “Técnico”: é mo- 
tivo de não menor alegria nem menor esperança no futuro da Instituição — porque na juventude elas 
residem, a alegria, a esperança, a força, a generosidade, a persistência, que fazem as Obras Grandes. 

Que seja grande a vossa Obra! 
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Palavras proferidas na Sessão de encerramento das comemorações dos 75 anos 
da Associação de Estudantes do I.S.T. por 


ABREU FARO 
Professor do I.ST. 


Julgo que Einstein nunca proferiu o lugar comum “'tudo é relativo”. Pelo menos no que respeita 
à Teoria da Relatividade. 

Pelo contrário. O Princípio da Relatividade liberta-nos de subjectivismo. Para qualquer observa- 
dor as leis da Física são as mesmas onde quer que esteja e como esteja. Até é possível escrevê-las sob 
a mesma forma. 

Assim, e sem que haja contradição, a Relatividade é coisa bem diferente do relativismo de um 
Protágoras e, mais nos nossos dias, de um mesmo relativismo de Kant. 

As Leis da Física, são da Física, mas tenho para mim que a verdade sendo verdade vai além do 
aspecto particular que reveste. Havendo humanização colhe-se luz das mais diversas fontes. 

Vem isto a propósito de quem pode apurar a verdade. 

Relaciona-se isto com o facto de quem e de que palavras nos podemos socorrer para testemunho. 

Fiéis ao Princípio da Relatividade qualquer um o pode fazer. 

No entanto, para que seja credível é necessário que ligue as palavras a juízos de realidade e de 
dignidade. E necessário que haja empenhamento, compromisso na afirmação de acções e resultados. 

Pela naturalidade e justeza que lhes são inerentes o enunciado das verdades primeiras corre sem- 
pre o risco de não se afirmar como conselho, de tal modo se confunde a aceitação dessas verdades 
com a prática que desse facto, e em coerência, deveria decorrer. E um dizer, um saber incontroverso 
de quase nenhum efeito. 

Isso, portanto, e infelizmente, não basta. 

Com esses enunciados universais se diz sim e se diz não e quase sempre o indesejável “tudo é 
relativo”, 

E necessário que as palavras do testemunho se enriqueçam com a legitimidade da experiência ou 
do recurso impecável aos primeiros princípios. 

Gozam destas propriedades as Notas Histórico-Pedagógicas sobre o Instituto Superior Técnico 
de Alfredo Bensaúde. 

“Nas notas que se seguem, exponho as ideias directrizes que nos guiaram, descrevo o que se pre- 
tendeu realizar e até que ponto foi possivel fazê-lo, comento os processos empregados no Instituto, 
em contraposição com os tradicionais nas nossas escolas, que pusemos de parte, e historio também 
algumas das dificuldades que foi preciso vencer para que a escola viitgasse... 

“Considero a redacção destas notas como o último acto da missão de que me encarregou o Go- 
verno Provisório..”” 

Que o observador pode ser qualquer prova-o o facto de muitas das passagens das Notas Histórico- 
-Pedagógicas poderem ter sido escritas por alunos. 

Omita-se o autor, saliente-se aqui e ali e não faltará quem diga que aquilo é obra de contestatário. 

As Notas de Bensaúde transcendem a realidade Técnico e são sem tempo. 

Há muitos anos emprestei a um amigo, mais idoso, culto mas que não conhecia o livro, o Conde 
de Abranhos... do Eça. | 

É espantoso como deixaram publicar! Embora escrito pelo Eça de Queiroz. É espantoso como 
é fácil ter-lhe acesso. 
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Assim, são essas notas de Bensaúde, plenas de sabedoria, de verdade, e em que se adivinham tem- 
pos difíceis de luta e sacrifício. 

Detesto histórias que à guisa de conclusão rematem com a moral da dita história, respeito a Moral 
mas contento-me com os factos, acho até antipedagógico que se saliente a moral. 

Leiam, leiam, as Notas Histórico-Pedagógicas de Bensaúde: alunos, professores. Todos aqueles 
que se preocupam com a Escola, com a Universidade, com o Ensino e Investigação. 

E difundam, difundam por todos os outros que sem grande preocupação e menor ocupação tam- 
bém falam da mesma Universidade. 

Em Bensaúde a vocação de professor era dominante, é um facto. 

Mas a força de que dispôs para realizar a sua obra; essa foi encontrá-la Bensaúde na crença que 
tinha na importância que reconhecia ao acto de ensinar, ao professor. 

E também na crença que, em abstracto, tinha nos alunos: 

“Há quem explique os fracos resultados obtidos em algumas das nossas escolas, como consequência 
da má qualidade dos alunos portugueses; é cómoda explicação ..:”. 

“A experiência de muitos anos permite-me asseverar que o estudante português é, em média, tão 
bom como o dos diversos países mais cultos da Europa..:” 

Com honra, prazer e alguma utilidade passei pelo Instituto de Alta Cultura. Assim foi-me possi- 
vel trazer até nós 9 estudantes japoneses que na Universidade de Sofia, em Tokio, seguiam a licencia- 
tura em Filologia Portuguesa. 

Foi agradável recebê-los. De facto foi inesquecível. 

Falaram um bom português. Japonês para mim seria impossível. Havia o inglês. Mas em portu- 
guês, em português era o ideal. 

Em vez de lhes falar de Portugal, perguntei-lhes pelo Japão. 

Julgo que é um país onde respeitam o professor? 

Entreolharam-se, uniram-se, falaram baixo e foi uma moça que me respondeu: 

Sim, é verdade: nos campos, nos teares, no teatro, na escola, todo aquele que ensina o que sabe 
é respeitado, é o professor. 

Ofereceram-me um livro, The Southern Barbarians, The First Europeans in Japan. 

Éramos nós. Lá chegámos nos meados do século XVI e ainda hoje nos agradecem a modernidade 
que então fomos capazes de lhes levar. Apropriaram-se do nosso modo de navegar, dos métodos de 
traçado das nossas cartas. 

O caso japonês é bem conhecido. 

Já antes. Mas sobretudo, depois da Primeira Grande Guerra, a partir de 1919, o Japão investiu 
poderosamente no ensino e na investigação. 

Visitou e enviou os seus filhos para centros de excelência do estrangeiro. 

Foi visitado, fez-se visitado por numerosos cientistas. 

Outra guerra, diversos alinhamentos políticos, até duas bombas atómicas lhes aconteceram. 

Mas as Leis da Física são sempre as mesmas e hoje aí estão. Aí está o Japão e o que faz. 

Ensinar é uma forma de dar. Se dar significa ficar com menos, nem sempre se dá. 

Alguém, desta casa, me perguntou se eu ensinava tudo quanto sabia. 

Tudo quanto souber em consciência e utilidade ensino, ensino tudo. 

Olhe que isso é perigoso! 

Lembrei-me, ocorreu-me então dizer que um bom professor assemelha-se, talvez, a um poço, dá 
de beber, esvazia-se um pouco, esvazia-se completamente, mas enche-se naturalmente, enche-se de novo 
e atinge o nível habitual no momento oportuno. 

Diferente é um depósito. Sem fontes próprias, que se enche custosamente e ainda mais custosa- 
mente se torna a encher. 

Mas pode haver outros motivos para não dar. Infelizmente há muitos. 
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Simples hipótese: uma ideia que sendo nossa se cede a outros. 

Imaginem Pasteur negociando, retardando a aplicação das suas investigações sobre a vacina. 

Talvez por isso, Pasteur dizia que não havia ciências puras e aplicadas, mas tão somente, aplica- 
ções da Ciência. 

E, sendo assim, como essa Ciência tem que ser pura, livre! Facultada a tudo e a todos a quem 
possa servir. 

Segredos de conveniência, competição nos lucros auferíveis, haverá sempre. Em todas as épocas, 
em todos os quadrantes. 

Mas cuidado! 

Deverá haver locais onde a generosidade, a isenção, a independência sejam o quotidiano e o pre- 
valente. Na Universidade em que me vejo deve ser assim. 

E para tanto não se define e ainda menos se quer uma Torre de Marfim. 

Na realidade a Universidade não pode alhear-se da sociedade onde se insere pela simples razão 
de, em princípio, a Universidade ser onde se sabe mais. 

Erro será não recorrer e não ajudar a Universidade. Grave será se, em qualquer medida, perder 
a sua independência. 

A Universidade ensina o que aprende através de uma actualização e investigação permanentes e 
continuadas. 

Mas o facto primeiro, aquele que deontologicamente lhe cabe, é ensinar. 

Ensinar âqueles que frequentam as licenciaturas e os seus cursos de pós-graduação. 

Ensinar a quem lhe pedir conselho. 

Colaborar na pesquisa de novas soluções. Sempre ensinar. 

No nosso entender o acto primeiro é ensinar e antes desse, anónima, difícil, desesperante por ve- 
zes, está a aprendizagem do que se vai ensinar, a investigação do facto novo com que se depara, o 
esclarecimento das contradições que se constroem. 

As leis de Newton são do domínio público, as equações de Maxwell não combram direitos de autor. 

O que torna legitimo o ensino destas matérias é tornar nosso o que foi descoberto e escrito por 
outros, é a assimilação séria e aprofundada do que se ensina. 


Mas também cuidado! 

Não confundamos essa ciência com a dispensa de se escreverem livros adequados. 

E estranho, um pouco bizarro, que se usem como livros de um curso, livros contemporâneos escri- 
tos por outros, nossos pares. 

Livro e professor devem constituir-se num todo harmonioso e em que cada parte é essencial. 

Não se confunde isso com bibliografia recomendada que deve ser abundante e acessível. 

E Bensaúde adverte: 


“Além do conhecimento do que ensinam os compêndios, o professor, para produzir bons discipu- 
los, precisa de outras aptidões, mais raras... 

« Pasteur, que era um fraco orador, se tivesse nascido português, teria passado a vida, provavel- 
mente, como figura secundária em algum dos nossos obscuros laboratórios. Num ambiente mais pro- 
pício, conseguiu ser um grande sábio e um grande professor, embora a retórica não fosse o seu forte. 
“Quando a sabedoria vai adiante a eloquência segue-a como fiel companheira” (Santo Agostinho)?” 


E sobre o aluno salienta: 

“O nosso aluno normal possui uma rapidez de compreensão igual e, às vezes, talvez superior à 
dos rapazes que frequentaram as escolas da Europa Central. E esta a maior das suas qualidades inatas; 
as que deveriam ser desenvolvidas pela educação faltam-lhe, porém, muitas vezes. O seu poder de assi- 
milação é quase sempre muito inferior ao dos estudantes cuja educação foi cuidada desde a infância. 
Pode dizer-se até que entre os nossos estudantes, a faculdade de assimilação costuma ser inversamente 
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proporcional à faculdade de compreensão; quanto mais facilmente compreendem, mais facilmente es- 
quecem. É que compreender e fixar o que se compreende são coisas diversas:” 

Sem Brito Camacho, Alfredo Bensaúde não teria fundado o Técnico, a Escola Nova em que o 
Instituto se converteu, 

Sem Duarte Pacheco, o Técnico nunca se teria desenvolvido como se desenvolveu. Isto, além do 
modo como Duarte Pacheco serviu a Engenharia Portuguesa e uma vez mais o Técnico. 

Sem condições mínimas, os espíritos fogem da carne frágil que os não agarra. 

Que se cuide do professor e da sua dignificação. 

Que se cuide do aluno e da sua dignificação. 

Deles depende o que for, mas outros, por dever de cargo, poderão e devem ajudar o que em mui- 
tos casos é de importância vital. 

A Associação dos Estudantes do Instituto Superior Técnico faz 75 anos. 

Criado por Decreto em 23 de Maio de 1911, o Técnico iniciou as aulas em 13 de Novembro. Seguiu-se 
a Associação dos Estudantes, fundada logo a seguir, 11 de Dezembro, exactamente como hoje. 

Desde sempre a Associação e o Instituto seguiram caminhos paralelos ou de cruzamento próximo 
com um destino comum: a perenidade do I.ST. 

Pela Associação passaram alunos distintíssimos que ocuparam lugares na Direcção, na Técnica, 
nas muitas e diversas áreas em que a associação se diversificou e aplicou. 

Passaram outros que não tendo sido alunos distintos foram, no entanto, hábeis gestores, impulsio- 
nadores e colaboradores da Associação. 

Foi sempre característica da Associação dos Estudantes do I.ST. o seu espírito fraterno de cama- 
radagem, entre si e com outras associações académicas. 

Foi sempre característica da Associação a sua inteira disponibilidade para com a Escola oferecen- 
do-se em generosa colaboração. 

De tal modo considero importante a acção das Associações de Estudantes que entendo por desne- 
cessário defender a sua existência. 

Além de tudo quanto realizou no campo pedagógico, cultural e desportivo. Além da realidade 
forte que é a TECNICA, além de tudo isso, a Associação foi e será sempre uma escola de dirigentes, 
não tenhamos dúvidas do facto, nem do que daí resultou e resultará. 

Houve momentos difíceis, de diálogo impossível e a maior parte das vezes impossibilitado. Mas 
houve sempre retorno amigo à colaboração. 

Política não é exactamente agitação nem outras coisas feias. Política, o Governo da Cidade, é uma 
coisa bela em que todos devem comparticipar para que não sejam governados como não querem. 

E uma coisa bela quando praticada nos exactos preceitos seguidos por uma Universidade indepen- 
dente e pura. 

As Associações foram e serão sempre origem natural de mentalidades robustas e críticas. De inter- 
venção possível e desejável nos sistemas educativos. 

Tenho a grande honra e acrescentado prazer de ser Director da Técnica. 

Na Técnica se publicaram, muitas vezes, as novidades científicas que, então, ainda incipientes, ga- 
nhariam vulto e extraordinária dimensão e influência. 

Essa consciência científica da Técnica deve-se aos professores do Instituto que nela colaboraram 
desde a primeira hora. 

A Técnica é uma revista portuguesa, de estudantes, escrita por professores e não exactamente de 
modo total para os estudantes mais novos. E um bom investimento, no entanto, para o futuro e uso 
desses mesmos estudantes. 

Numa tendência para certo derrotismo, numa vaidade sem ter de quê, nem sempre analisada até 
às raizes onde se alimenta, nem sempre se acorre à Técnica. É pena. É mau para o pensamento cienti- 
fico português. É mau exemplo que não ajuda a frutificar outras revistas similares. 
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No entanto, a Técnica viaja, pedem-nos separatas de além-fronteiras. 

Termino. 

Associação é uma palavra bela, plena de intenções. 

Estudante, ser estudante, é belo e sem idade. 

Os jovens terão sempre por si o futuro que postulam e invadem e para o que, em boa verdade, 
não há experiência que invalide os sonhos com que, no presente, o controem. 

A vida continua-se com a chegada dos que nascem. 

Estudante na escola é episódico e efémero. 

Que a Escola e a sua Associação consigam formar estudantes de sempre e associados de sempre. 

Isso significa a perenidade da trilogia: 

Instituto Superior Técnico. 

Associação dos Estudantes do I.S.T. 

Associação dos Antigos Alunos do I.ST. 

Tão só isso! 
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Palavras proferidas na Sessão de encerramento das comemorações dos 75 anos 
da Associação de Estudantes do I.ST. por 


FRANCISCO SANCHES 
da Direcção da AEIST 


Celebrar 75 anos de vida de uma organização, é necessariamente, reflectir no passado, analisar 
o presente e projectar o futuro. | 

Quando se trata duma Associação fundada e gerida por jovens universitários, junta-se ainda a 
obrigação de realizar cada vez mais e melhor. 

Em 11 de Dezembro de 1911 os objectivos estatutários da Associação de Estudantes do Técnico 
eram os seguintes: 


“a) Difundir a instrução entre os associados, promovendo excursões e visitas de estudo, conferên- 
cias literárias e científicas e especialmente as que se relacionam com o estado financeiro, indus- 
trial, comercial e agrícola do nosso país e tudo o mais que possa concorrer para o aperfeiçoa- 
mento moral e físico dos seus sócios. 

b) Procurar obter regalias que tornem as condições dos associados mais fáceis.” 


Como se nota são objectivos que transportados para os nossos dias não perderam actualidade. 

Pensamos que ao longo do tempo a Associação de Estudantes do IST tem sabido cumprir as ideias 
de solidariedade e de ajuda aos estudantes do Técnico. 

Depois de uma fase inicial em que acompanhou uma certa indefinição da escola, em termos de 
fins e instalações, a Associação de Estudantes passou a dispor de pavilhão próprio, passo importante 
para a sua consolidação e expansão. 

Seguiram-se no tempo fases opostas que ““adormeceram”” um pouco a acção da Associação. 
Referimo-nos ao controlo dos corpos gerentes exercido pelo Estado e ao fervor revolucionário após 
Abril 1974. Foram fases das quais é importante tirar os ensinamentos necessários para não cair nos 
mesmos erros. 

Actualmente, a Associação de Estudantes do Técnico apresenta, uma estrutura dividida por várias 
áreas de intervenção: Associativa, Federativa, Desportiva e Pedagógica subdividida por secções ama- 
doras que promovem e sustentam as actividades extracurricular dos alunos do I.ST. Efectivamente, 
é no parque desportivo e nas suas restantes instalações que é dinamizada a actividade recreativa/cultu- 
ral e desportiva dos estudantes. Infelizmente, as cargas horárias continuam a não permitir um acesso 
regular às actividades extracurriculares. 

Se no plano interno se reconhece que é necessário melhorar os serviços prestados, é no plano ex- 
terno que se consolida a posição da maior e mais importante Associação de Estudantes da Academia 
de Lisboa. Não é por acaso que as duas estruturas federativas mais recentes, Associação Académica 
de Lisboa e a Rádio Universidade Tejo são lideradas por representantes da nossa Associação. 

Outro vector importante de actuação é o da intervenção pedagógica contribuindo para a melhoria 
do ensino no I.S.T. Pensamos que será aqui que a Associação deverá intervir no futuro ainda com mais 
empenho e maior participação de todos os estudantes. 

Há que ultrapassar a velha luta aluno-professor avançando para a resolução dos problemas atra- 
vés de ideias e projectos que indiquem que a nossa função não é a da contestação pura e simples. Apostar 
num ensino melhor é a meta que a Associação de Estudantes do I.S.T. deve atingir o mais depressa 
possível. No entanto, isto só se consegue se houver a participação de todos, através de sugestões e críticas. 
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A Associação de Estudantes do Técnico será aquilo que todos os Estudantes desejarem, pois, há 
sempre a possibilidade de mudar os seus dirigentes e avançar com alternativas. 

Continuar a obra iniciada em 11 de Dezembro de 1911 não é tarefa fácil. Mas se todos quiserem 
teremos uma Associação melhor e mais participada. 
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Palavras proferidas na Sessão de encerramento das comemorações dos 75 anos 
da Associação de Estudantes do I.ST. por 


VENÂNCIO ROSA 
da Comissão Organizadora das Comemorações 


O motivo que nos trouxe aqui tem algo de especial, pois a Associação de Estudantes do I.S.T. 
faz hoje 75 anos. 

Desde 11 de Dezembro de 1911 — data da sua fundação — até aos nossos dias esta Associação 
de Estudantes tem sido um marco na vida académica do País, nunca deixando que a sua voz esmoreça 
ou seja abafada sempre que estão em causa os direitos que a todos os estudantes devem assistir, e espe- 
cialmente na cidade de Lisboa onde, pelo seu dinamismo muito tem contribuído para estreitar os laços 
que unem as diversas faculdades onde temos, como exemplo recente, o facto de ser uma das Associa- 
ções de Estudantes que impulsionou a formação da A.AL.L. 

Por tudo isto e porque são poucos os organismos deste tipo que se conseguem manter durante 
tantos anos com tão forte personalidade, não quisemos deixar de comemorar convenientemente estas 
Bodas de Diamante. 

Estas comemorações tiveram como principais objectivos: 


— mostrar O que os alunos fazem, hoje, em cada curso, no Técnico e em que termos a indústria 
e o País podem contar com os futuros engenheiros; 

— abrir as portas para o exterior, especialmente para a cidade; 

— promover o convívio entre universitários das diversas faculdades que, como todos sabemos, é 
bastante deficiente, devido a estas estarem muito dispersas — facto que muito lamentamos. 


Para conseguirmos estes objectivos realizámos um vasto leque de actividades que tiveram início 
no dia 17 de Novembro e se prolongaram até hoje, das quais passo a fazer uma breve retrospectiva. 

De 17 a 22 de Novembro decorreu a exposição “A Engenharia no I.S.T.”” que, contando com equi- 
pamento de laboratório dos vários cursos e com muitos trabalhos de final de curso executados pelos 
alunos mostrou a todos os visitantes, que foram muitos, as grandes potencialidades que o engenheiro 
licenciado no Técnico tem, tanto para a indústria como para a investigação. 

De 22 a 30 de Novembro realizaram-se diversas actividades das quais passo a referir algumas: 


— Torneios inter-faculdades de basquetebol, voleibol, futebol-5, andebol, pólo aquático, uma prova 
de atletismo e mesmo um torneio de xadrez, o que movimentou centenas de universitários. 

— Cinema — já é tradição passarem-se bons ciclos de cinema no I.S.T. Esta foi mais uma ocasião 
para a confirmação deste facto, contando-se com três filmes que tiveram uma afluência de pú- 
blico além das espectativas. 

— Teatro — não sendo muito comum o gosto pelo teatro no I.ST. o professor Listopad e o Grupo 
de Teatro da UÍT.L. vieram mostrar-nos como se faz teatro universitário. 

— Noite de Coimbra — é quase um hábito contarmos nas nossas festas e comemorações com 
os colegas de Coimbra cuja tradição marca a vida académica do País. Contámos nessa noite 
com a presença da Tuna, Orfeon e um grupo de fados e cantares da Associação Académica 
de Coimbra que nos deliciaram com a boa música e cantares de Coimbra e contou com um 
público caloroso que esgotou a lotação da sala. 

Melhor ainda foi o Convívio que se estabeleceu entre todos os elementos, cerca de 100 e os 
alunos do Técnico que com eles contactaram. 
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— Baile do Engenheiro — com este baile quisemos reatar a tradição dos bailes de gala nesta casa, 
mas teve como objectivo principal o estreitar das relações entre estudantes, professores e anti- 
gos alunos. 


No dia 6 de Dezembro realizámos um espectáculo na Aula Magna, pois, não quisemos deixar pas- 
sar estas comemorações sem fazer referência à música popular portuguesa, por isso, organizámos um 
espectáculo com o cantor Fausto que cantou com a sala cheia, verificámos com alegria, que a maioria 
eram estudantes. 

A Comissão Organizadora destas comemorações e a D.A.E.I.ST. querem agradecer a todas as pessoas 
que ajudaram na organização e colaboraram connosco e, muito especialmente, aos membros da Comis- 
são de Honra, que sem o seu apoio e avalo não tinha sido possível a realização de todas estas actividades. 

Queremos também aproveitar esta oportunidade para prestar homenagem a todos os dirigentes 
associativos que ao longo destes 75 anos trabalharam para que a nossa Associação de Estudantes seja 
o que é hoje e, principalmente, aqueles que não querendo deixar que o lema destes 75 anos “Construir 
um futuro” se ficasse pela sua vida académica, se distinguiriam, por isso, na vida pública nacional. 
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FÍSICA DE ALTAS ENERGIAS: Actividade recente 
do Grupo de Física Experimental de Altas Energias 


de Lisboa 


José Mariano GAGO*, Amélia MAIO**, 
Mário PIMENTA, João VARELA* 


RESUMO 


Apresenta-se sucintamente algumas das experiências 
de física de altas energias no CERN nas quais participa 
o Grupo de Física Experimental de Altas Energias de Lis- 
boa, indicando, para cada uma delas, os objectivos pro- 
postos, o equipamento utilizado, e os resultados alcan- 
çados. 

Apresenta-se, igualmente, um trabalho de desenvol- 
vimento tecnológico recente no domínio das fibras ópti- 
cas cintilantes, motivado pela necessidade de novos de- 
tectores em física de altas energias. 


1. INTRODUÇÃO 


O Grupo de Física Experimental de Altas Energias 
de Lisboa é um colectivo de trabalho formado por in- 
vestigadores de várias instituições para actividades de pes- 
quisa no domínio da física das partículas elementares. O 
Grupo realiza parte do seu trabalho em Lisboa, nomea- 
damente a concepção de experiências, desenvolvimento 
tecnológico, desenvolvimento de “software” e andlise de 
resultados, executando no CERN o que diz respeito à pre- 
paração directa das experiências e à tomada de dados, 
em colaboração com grupos de outros países'”. 


* Instituto Superior Técnico (Departamento de Física) e CEMC (INIO). 
** Faculdade de Ciências de Lisboa (Departamento de Física) e CFNUL 
(INTO). 


“ A composição do Grupo de Física Experimental de Altas Energias de 
Lisboa é a seguinte (Novembro 1985): a) Investigadores: Conceição Abreu 
(FCL), Gaspar Barreira (CFNUL), Paula Bordalo (IST), Augusto Casaca (IST- 
-INESC), José Mariano Gago (IST), Amélia Maio (FCL), Mário Pimenta 
(IST), Sérgio Ramos (IST), João Varela (IST); b) Estagiários: Carlos Cruz 
(IST), Paulo Filipe (IST), Luis Peralta (FCL), Pedro Vaz (IST). O GFEAEL 
é coordenado por José Mariano Gago e funciona em Lisboa, no Complexo 
H do INIC (Instituto de Física e Matemática), Av. Prof. Gama Pinto, 2, 1699 
Lisboa Codex, Tel.: 773325/773338/774297, Telex 62593 HFM P. 


Original recebido para publicação em 20/11/85. 


Este artigo visa apresentar algumas das actividades 
recentes do Grupo. 

Falaremos assim de três grandes tópicos. Em primeiro 
lugar, focaremos uma experiência em preparação, cujo pe- 
ríodo de tomada de dados terá lugar em finais de 1986 
e início de 1987, destinada à investigação da possível ener- 
gia. Esta experiência, denominada NA38, detectará pa- 
res de muões emitidos nessas colisões e medirá ainda a 
energia neutra emitida. Sucede à experiência NAIO, na 
qual também participam físicos portugueses, onde foram 
feitas descobertas importantes sobre o comportamento dos 
quarks e gluões no interior da matéria hadrónica. 

Em segundo lugar, descreveremos uma experiência já 
realizada e concluída, destinada a esclarecer a natureza 
do pico na secção eficaz de retrodifusão das colisões 
antiprotão-protão, na região dos 180º, aparentemente idên- 
tico ao crescimento que se observa em óptica nas secções 
eficazes de retrodifusão da luz pelas gotas de dgua e que 
é responsável pelo efeito chamado de “auréola”. A ex- 
periência permitiu mostrar que esse pico da secção efi- 
caz antiprotão-protão não é de facto explicável como 
efeito de auréola de tipo óptico; permitiu ainda chamar 
a atenção para um comportamento em energia das sec- 
ções eficazes de certos canais, cujo esclarecimento exigirá 
outros dados experimentais. 

A mesma experiência WA 74 conduziu a um resultado 
inesperado em domínio diverso daquele para que fora ini- 
cialmente projectada. Com efeito, foi possível estudar, com 
precisão nunca atingida anteriormente, a polarização dos 
lambdas produzidos na região extrema para a frente (x 1) 
com feixe incidente de K”. A polarização medida (que 
é em última análise, a do quarks difundido a pequeno 
ângulo na matéria hadrónica) cresce fortemente com o 
momento transverso do lambda produzido. A razão deste 
importante fenómeno (só agora descoberto) não é co- 
nhecido. 

Em terceiro e último lugar, descrevemos o trabalho 
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desenvolvido no domínio da utilização e do estudo de fi- 
bras plásticas cintilantes, cujo sucesso parece assegurado 
para a realização de certo tipo de detectores, nomeada- 
mente de calorímetros electromagnéticos e de hodoscó- 
pios de cintilação de grande resolução. 


2. A EXPERIÊNCIA NA38 (em preparação): 
PROCURA DE PLASMA DE QUARKS 
E GLUÔESO 


A experiência NAIO sondou a estrutura da matéria 
hadrónica estudando as colisões de piões de 200 GeV so- 
bre um alvo de tungsténio onde são produzidos pares 
de muões de carga contrária. Estes são coerentemente 
interpretados como o resultado da aniquilação de um 
quark com um antiquark, respectivamente de um nucleão 
no alvo e de um pião incidente, num fotão virtual mas- 
sivo que se materializa num par de muões. Os resulta- 
dos confirmam que os hadrões são partículas com es- 
trutura, dão a distribuição de impulsão dos quarks no 
seio dos piões e nucleões e pôem em evidência alguns 
dos aspectos da dinâmica da interacção forte. 

Esta visão da matéria composta de quarks estruturou- 
-se de facto com base nos resultados de numerosas ex- 
periências diferentes, embora nenhuma das tentativas rea- 
lizadas para observar os quarks no estado livre tivesse 
resultado. Actualmente julga-se compreender este facto 
como o resultado da natureza da interacção forte que 
obrigaria os quarks a manterem-se confinados no inte- 
rior dos hadrões. 

A matéria pode, no entanto, organizar-se diferente- 
mente quando submetida a condições extremas de tem- 
peraturas e pressão, tais como as que existiram nos pri- 
meiros instantes da criação do Universo (A Grande 
Explosão) e as que existem actualmente no interior das 
estrelas de neutrões. De facto, a teoria que os físicos dis- 
põem para descrever a interacção entre os hadrões (Cro- 
modinâmica Quântica), quando extrapolada para den- 
sidades de energia superiores à observada na matéria 
nuclear (e, = 0.15 GeV/fm”) prevê a existência de uma 
transição de fase para um outro estado da matéria, o 
plasma de quarks e guões desconfinados. Quando a den- 
sidade de energia atinge valores elevados, quer devido 
a uma compressão extrema, caso das estrelas de neutrões, 
quer por via de grandes energias de agitação térmica 
(temperatura), caso dos instantes iniciais do Universo, 
o campo de cor (carga da interacção forte) sofre um 
efeito de blindagem ““tipo Debye”, de forma que os ha- 
drões perdem a sua identidade e os quarks e guões co- 


* A colaboração cientifica responsável pela experiência NA38 é a seguinte: 
Bergen, CERN, Clermont-Ferrand, Ecole Polvthechnique, Lisboa (INIC), 


Ivon, Neuchâtel, Orsay, Strasbourg, Valência. 


Os signatários portugueses desta proposta de experiência são: €. Abreu, 
G. Barreira, P. Bordalo, A. Casaca, J. M. Gago, A. Maio, M. Pimenta, S. Ra- 


mos e J. Varela. 


loridos passam a formar um gás de partículas quase- 
livres. As estimativas existentes indicam que a transição 
de fase se dá para uma densidade de energia da ordem 
de 2 GeV/fm”, o que no plano temperatura-densidade de 
número bariónico corresponde à linha esboçada na 
fig. 1.1. 


Ech 


Em o 
2 2 


a. 


Fig. 11.1 a) Representação esquemática de dois núcleos 
observados no referencial do Centro de Massa antes da 
colisão. As dimensões longitudinais são da ordem de 
1 fm; b) Região de formação de matéria após a colisão. 


Num futuro próximo abre-se a possibilidade de re- 
produzir em laboratório as condições iniciais do Uni- 
verso, embora, como é óbvio à escala das dimensões nu- 
cleares. De facto, preparam-se actualmente, no CERN 
experiências destinadas ao estudo das colisões com iões 
de oxigénio acelerados a 200 GeV/nucleão no SPS (Su- 
per Proton Syncrotron). Nas colisões mais violentas, 
estima-se que sejam produzidas da ordem de 400 parti- 
culas e que se atinjam densidades de energia da ordem 
de 2 GeV/fmº. Para compreender o que se passa quando 
dois nucleões colidem a alta energia é necessário fazer 
apelo à teoria da relatividade. Observados no referen- 
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Fig. 1.2 Representação da frase hadrónica e do plasma 
de quarks e gluões no diagrama temperatura-densidade 
de número bariónico. 
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Fig. II.3 Espectrómetro de muões da experiência NA.10. 
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Fig. 11.4 Detectores na região do alvo. 
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cial do centro de massa da interacção, os dois núcleos, 
por contracção relativista, apresentam a forma de duas 
bolachas movendo-se em direcções opostas (fig. 11.2). As 
dimensões longitudinais dos núcleos, que no limite da 
velocidade da luz se anulariam são limitados a 1 fm pelo 
princípio de incerteza. O tempo de colisão, é portanto, 
da ordem de 1 fm/c. Nas colisões são criadas partículas 
num tempo próprio da ordem das suas dimensões, 
| fm/c. No referencial do centro de massa este tempo é 
dilatado relativisticamente resultando, portanto, consi- 
deravelmente superior ao tempo de colisão. Desta forma 
os núcleos são transparentes no momento de interacção, 
ocorrendo a criação de matéria (quarks e gluões) em ins- 
tantes posteriores, na região delimitada pelos núcleos em 
afastamento. Se a densidade de energia nesta região for 
inferior à densidade crítica, os quanta de matéria 
organizam-se en hadrões constituindo um gás de parti- 
culas que se expandirá até ser detectado pelo equipa- 
mento experimental. No caso contrário não há lugar para 
a imediata condensação hadrónica, de forma que os 
quarks e gluões permanecem livres formando um gás su- 
jeito a dois processos em competição, um de termaliza- 
ção conducente à uniformização da temperatura, e ou- 
tro de expansão, essencialmente longitudinal, governado 
pelas equações da hidrodinâmica relativista. Em tempos 
posteriores do processo de expansão o gás arrefece, a den- 
sidade de energia desde abaixo de ponto crítico, e os 
quarks recombinam-se dando origem aos produtos fi- 
nais da reacção. 

Apercebemo-nos facilmente de que a complexidade 
do processo de colisão de dois núcleos relativistas obs- 
curece a eventual produção do plasma de quarks e gluões, 
tornando delicada a sua identificação e caracterização. 
Com efeito, a grande maioria dos produtos da reacção 
reflecte, essencialmente, a física que ocorre na frente de 
hadronização ao longo da expansão do plasma. No en- 
tanto, o gás de quarks e antiquarks deve ser a sede de 
abundante produção de fotões, distribuídos segundo a 
lei de Planck se aquele se encontrar termalizado. Estes 
fotões, tanto reais como virtuais, são um verdadeiro ter- 
mómetro indicador das características do plasma e da 
natureza da transição de fase. 

Do ponto de vista experimental, a detecção de fo- 
tões reais provenientes do plasma é seriamente dificul- 
tada pela presença dos fotões de desintegração dos piões 
neutros, abundantemente produzidos no processo de har- 
monização. Em contrapartida, os fotões virtuais mate- 
rializam-se em pares de muões de fácil detecção e me- 
dida. A experiência NA38 no CERN visa, assim, a 
medida dos pares de muões de massa superior a 0.5 
GeV/c? produzidos nas colisões oxigénio-urânio, acom- 
panhada da medida da distribuição de energia hadró- 
nica por forma a determinar o grau de inelasticidade das 
colisões. 

A detecção de pares de muões é feita pelo espectró- 
metro da experiência de NAIJO (fig. 11.3). O princípio de 
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detecção baseia-se na diferença das probabilidades de in- 
teracção com a matéria dos muões e dos hadrões. Um 
bloco de carbono de 4,8m após o alvo absorve a totali- 
dade dos hadrões produzidos conjuntamente com os 
muões, que são desta forma isolados e medidas pelos 
detectores colocados a jusante. A medida da impulsão 
é feita por deflecção das trajectórias num campo mag- 
nético toroidal, as quais são reconstruídas a partir da 
informação recolhida em câmaras de fios. Os aconteci- 
mentos interessantes são sinalizados por lógica electró- 
nica de grande rapidez, com base nos sinais fornecidos 
por um conjunto de 800 cintiladores reunidos em 4 ho- 
doscópios. 

A medida dos hadrões é realizada por um detector 
compacto, instalado entre o alvo e o filtro de muões, es- 
pecialmente concebido para esta experiência (fig. 11.4). 
Um calorímetro electromagnético mede a distribuição de 
energia dos fotões produzidos, os quais resultando es- 
sencialmente do decaimento 7º>y +y nos dão uma 
medida dos piões neutros. Os 36 módulos do caloríime- 
tro são constituídos por camadas sucessivas de chumbo 
e fibras ópticas cintilantes de 0.5 mm de espessura. A luz 
recolhida nas fibras cintilantes é proporcional à energia 
dos fotões incidentes iniciadores de cascatas electromag- 
néticas no interior do calorímetro. Os sinais luminosos 
são conduzidos por guias de luz até fotomultiplicado- 
res onde são convertidos em impulsos eléctricos, poste- 
riormente digitalizados por conversores analógico-di- 
gitais. 

O alvo de urânio subdividido em 20 segmentos de 
| mm de espessura espaçados de 1 cm, é envolvido por 
25 cintiladores anelares de 0.7 cm de diâmetro. A altura 
de impulso medida nos cintiladores é proporcional ao 
número de partículas carregadas que os atravessam. Este 
sistema cumpre duas funções, por um lado mede a dis- 
tribuição angular das partículas carregadas e por outro 
identifica eventuais reinteracções de fragmentos nucleares. 

O feixe incidente, cuja intensidade é de 10º iões/seg, 
contém uma pequena proporção (1%) de iões de cálcio, 
para além de fragmentos de menor número atómico (Be, 
alfa...) É, portanto, necessário dispôr de um sistema ca- 
paz de medir intensidades elevadas e de identificar os 
iões incidentes. Estas funções são realizadas complemen- 
tarmente por um hodoscópio de 28 cintiladores colocado 
no percurso do feixe e por dois detectores de quartzo 
colocados à entrada e à saída do alvo. 

O conjunto dos sinais analógicos fornecido pelo equi- 
pamento adicional (calorímetro, detector de alvo, detector 
de feixe) é analisado em amplitude por electrónica apro- 
priada e transferido, conjuntamente com os sinais lógi- 
cos provenientes do detector de muões, para um sistema 
de computadores Nord 100/500, onde são registados em 
banda magnética para análise posterior. 

O Grupo de Física Experimental de Altas Energias 
de Lisboa, (GFEAEL) participou activamente na elabo- 
ração da proposta de experiência NA38, recentemente 
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aprovada pelo Comité de Pesquisa do CERN, nomea- 
damente os estudos de adaptação do espectrómetro de 
muões às novas condições experimentais. Durante o pró- 
ximo ano, a GFEAEL é responsável pelo desenvolvi- 
mento dos programas de aquisição de dados em linha, 
bem como dos programas de simulação e por testes do 
calorímetro electromagnético. O GFEAEL é também res- 
ponsável, em colaboração com o INESC, pela constru- 
ção de uma interface destinada a compatibilizar os dois 
sistemas de aquisição utilizados (RMH-ROMULUS). 


3. A EXPERIÊNCIA WA74 (concluída) 


A experiência WA74 realizou-se entre Fevereiro e 
Maio de 1982 no espectrómetro Omega na zona oeste 
do CERN, urilizando um feixe não separado (e”, m7, k- 
p”) e um alvo de hidrogénio líquido (fig. II.1.). 

O espectrómetro Omega dispõe de um iman supra- 
condutor com um campo magnético vertical ( - 1.8T) que 
permite boas aceitâncias geométricas. 

As trajectórias das partículas carregadas foram re- 
construídas, na zoná para a frente, por um sistema de 
quinze câmaras proporcionais e duas câmaras de deriva 
e, na região de grandes ângulos por dezasseis câmaras 
proporcionais colocadas paralelamente ao alvo. 

As partículas rápidas difundidas para a frente atra- 
vessam um ou dois Cerenkovs de limiar, que permitem 
estabelecer, para determinados intervalos em momento, 
hipóteses de massas, e dois ou três hodoscópios verticais. 

O trigger da experiência exigia a existência de uma 
partícula rápida positiva de impulsão próxima da dos 
feixes e foi elaborado com a ajuda da informação reco- 
lhida por um sistema de leitura rápida instalada num 
plano duma das câmaras proporcionais. 

Este trigger permitiu seleccionar candidatos às reac- 
ções a dois corpos 7” p>p7r”, pp>7m* 77, k-p>pk- e 
pp pp e ainda estudar produção inclusiva de lambdas. 

Nos parágrafos seguintes apresentam-se os principais 
resultados obtidos nestes dois tipos de reacções: 


a) Estudos de reacções a dois corpos 


Na reacção m- p—p7r” a secção eficaz diferencial 
do/du foi determinada entre O<u<1 GeV?/c?. Os valo- 
res encontrados nesta recção enquadram-se nos valores 
publicados para momentos incidentes entre 6 e 90 GeV/c 
que são globalmente explicados num modelo de Regge 
pela troca de uma única trajectória com 
a = —(0.06 + 0.05) + (0.50 + 0.26)u. 

O valor encontrado para a secção eficaz na reacção 
pp>m*am” a 12 GeV/c permitiu, comparando com re- 
sultados publicados a 5 e 6GeV/c, estabelecer uma de- 


* A colaboração responsável pela experiência WA74 é a seguinte: CERN, 
Universidade de Paris VI-LPNHE, Collêge de France, Lisboa (INIC). 
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Fig. HI.1 Dispositivo Experimental experiência WA 74. 
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Fig. HI2 Secção eficaz diferencial na região para trás 
(u=O) em função da energia no centro de massa na 
reacção pp—pp. 
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Fig. 11.3 Secção eficaz diferencial dO/dt para vários 
momentos incidentges na reacção pp>pp. 


pendência em energia análoga, dentro dos erros, à veri- 
ficada na sua reacção cruzada (mr ppm”). Este 
resultado comprova pela primeira vez, para estas reac- 
ções e numa gama de energias significativa, uma previ- 
são geral dos modelos de Regge pela qual as trajectó- 
rias trocadas em reacções cruzadas são idênticas. 

Os limites superiores encontrados a 8e 12 GeV/c na 
secção eficaz para trás na reacção pp—>pp indicam a 
manutenção do regime de forte dependência em energia 
estabelecido anteriormente nesta reacção para momen- 
tos incidentes inferiores a 6.2 GeV/c (fig. 111.2). De igual 
modo estes limites são, quando comparados com os re- 
sultados encontrados na mesma reacção e energias na 
região de grandes ângulos, compatíveis com o cresci- 
mento da secção eficaz diferencial do/du para u= 0 
GeV/c? como tinha sido verificado a 5 e 6.2 GeV/c 
(fig. 111.3). 

Na reacção k-p—pk- os valores de sensitividade al- 
cançados com o feixe k não permitem que os limites 
superiores obtidos na secção eficaz para trás tenham sig- 
nificado apreciável. 

A coexistência, em reacções exóticas, de uma forte 
dependência em energia da secção eficaz para trás com 
o crescimento, na mesma região, da secção eficaz dife- 
rencial do/du para u = O Gev/c” corresponde a uma mu- 
dança de regime na dependência em energia na secção 
eficaz para trás passando-se de uma forte dependência 
para uma dependência suave. Este comportamento é as- 
sociado, em modelos de Regge, ao aparecimento de uma 
trajectória dominante. Note-se, porém que, mesmo para 
os canais não exóticos, não existe qualquer modelo teó- 
rico que explique a zona de forte dependência em ener- 


gia. Esta dependência é aliás análoga à encontrada nos 
canais exóticos na região para trás bem como nas reac- 
ções exclusivas na região de grandes ângulos. 

O estudo exploratório do comportamento das sec- 
ções eficazes na região para trás em reacções exóticas 
a energias mais elevadas foi efectuado nesta experiên- 
cia. O problema físico subjacente, compreensão dos me- 
canismos dominantes nas zonas de forte dependência em 
energia, parece-nos interessante e merecedor dum estudo 
sistemático. Para tanto é indispensável um feixe separado 
e intenso de p e/ou k- com um momento incidente en- 
tre 10 e 20 GeV/c que permita sensitividades de - 1 acon- 
tecimento por picobarn. De momento, não se encontram 
disponíveis feixes com estas características. Contudo, a 
existência de alguns projectos em estado de elaboração 
muito diferente (nova geração de acumuladores p, 
TRIUME, HIPS, SUPER-LEAR...), permitem ter espe- 
rança que num horizonte de dez anos se possa modifi- 
car significativamente este panorama. 


b) Polarização dos lambdas emitidos na região de frag- 
mentação do K” 


Nesta componente da experiência WA74 a identifi- 
cação da partícula do feixe foi obtida usando dois de- 
tectores Cerenkov na linha de feixe: um detector de li- 
miar e um detector diferencial (CEDAR) ajustado à 
identificação de K”. O “trigger” que selecciona um pro- 
tão rápido emergente da colisão aceita lambdas no modo 
de desintegração A-p7”. A aceitância dessa selecção 
cresce de 10% a x = 0.6 até 70% a x = 1.0, onde x é de- 
finido como x=Pt%Pf |, razão entre o momento lon- 
gintudinal do lambda no centro de massa da colisão K” 
p e o momento longintudinal máximo possível nesse re- 
ferencial. 


Após reconstrução e selecção, foram obtidos 16 630 
lambdas a 12 GeV/c, mais de 90% dos quais com x>0.6. 
Apenas estes foram, finalmente, utilizados na análise. 

A polarização dos lambdas pode ser estimada por 
<aP>=<cos0>/<cos'0>, onde «=0.642 e cos= 
=(K xA).Pm/ |(K x A)|.| Pam |. Nesta fórmula, 
Pema é O momento do protão de desintegração no refe- 
rencial próprio do lambda. 

Após uma análise exaustiva destinada a corrigir e mi- 
nizar todos os possíveis erros sistemáticos introduzidos 
pela aparelhagem e pelos critérios e processos de recons- 
trução e selecção, obtiveram-se, no que diz respeito à po- 
larização dos lambdas, os resultados que se mostram nas 
figuras 11.4 e JII.5. 


Resumidamente, esses resultados mostram que os 
lambdas produzidos na região de fragmentação do K' 
são polarizados na direcção K x A. Essa polarização 
aumenta com o momento transverso do lambda (p,) e 
é insensível a x na região considerada. A dependência 
da polarização com p, pode ser ajustada à fórmula 
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Fig. HI.4 Polarização do 
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Fig. III.5 Polarização do 
Lambda em função do 
momento transverso do 
Lambda a 12 e 16 GeV de 
momento incidente. 
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P (p) = Bp, com o resultado B = 0.66+0.03 (x” = 6.0/ 
/ND = 3). 

A importância deste resultado, que é inteiramente 
novo, reside no facto de a polarização dos lambdas na 
região de fragmentação do K'” ser, em última análise, a 
polarização adquirida pelo quark s, inicialmente não po- 
larizado, ao difundir-se (e recombinar-se) na matéria ha- 
drónica. 

Um processo de polarização por difusão é conhecida, 
por exemplo, em electrodinâmica quântica, aquando da 
difusão de electrões não polarizados por um potencial 
coulombiano ou por um núcleo; a polarização é um re- 
sultado, nesse caso, dos diagramas de 2º ordem de QED 
(polarização de Mott). Podia, pois, legitimamente, per- 
guntar-se se comportamento idêntico se não previria em 
cromodinâmica quântica. A resposta é negativa, ao me- 
nos no que diz respeito aos diagramas de 2º ordem: a 
dependência angular da polarização e a dependência em 
energia, nesse caso, são incompatíveis com os nossos re- 
sultados experimentais que permanecem, assim, e de mo- 
mento, sem explicação. 


4. FIBRAS ÓPTICAS CINTILANTES DE PLÁSTICO 
(desenvolvimento técnico para detectores de partículas) 


As fibras ópticas cintilantes de plástico têm sido de- 
senvolvidas para detecção de radiação electromagnética 
e em particular para calometria electromagnética em Fi- 
sica das Altas Energias. Por essa razão, o estudo das suas 
propriedades ópticas e de cintilação, e a sua resistência 
à radiação (U.V., y, n, e”) é da maior importância. 


As fibras ópticas cintilantes'” (f.o.e.p.) em estudo, 
são constituídas por um corpo central cilíndrico de po- 
liestireno revestido de uma fina camada de acetato de 
polivinilo. O poliestireno é dopado com 1% de Butil- 
-PBD, a componente cintilante, e com 10-* de Demetil- 
-Popop, o 2º soluto, cuja função é deslocar a emissão 
no UV. próximo para o visível. O D-Popop absorve a 
luz de comprimento de onda À - 375 mm e reedite no azul 
de c.d.o. À = 430 mm (fibra azul). Como o 2º soluto or- 
gânico outros compostos estão sendo ensaiados, como 


para o foto - 
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Fig. IV.1 Esquema da produção de luz e sua transmis- 


são ao longo da fibra. 


n,=1,59; n, = 1,46; AN = 0,62. 


sejam o BBOT, o BDB e o K27. Em particular o K27 
que emite no verde (fibra verde) têm-se apresentado como 
possível alternativa ao D-Popop. 

O esquema de produção de luz e sua transmissão de- 
vida ao fenómeno óptico de reflexão total, estão repre- 
sentadas na fig. IV.l. 


Como se indica na fig. 1, n, é o índice de refracção 
do corpo central e n, é o índice de refracção do reves- 
timento. Para estas fibras os valores dos índices de re- 
fracção são n,= 1.59 en,= 1.46 e o valor da abertura 
numérica é de AN = 0.62. No momento actual são fa- 
bricados fibras com diâmetros que vão de 0.5 mm a 2mm 
e a espessura do revestimento correspondente varia de 
Sum a 20 um. 

O objectivo deste projecto é pois, o estudo 


1) das propriedades de cintilações das fibras ópticas 
de plástico cintilante, 1.e., o rendimento na pro- 
dução de luz quando a fibra é atravessada por ra- 
diação electromagnética, e 

1) das características de transmissão dessa luz atra- 
vés da fibra, medindo o comprimento de atenua- 
ção L, da fibra, i.e., o comprimento de fibra 


* Desenvolvidas em DP,PE — STIPE, SACLAY, FRANÇA. Em Lisboa 
o trabalho aqui referido foi realizado por Carlos Cruz (CFMC), Amélia Maio 
(CFNUL) e Luis Peralta (CFNUL), sob a direcção da Prof? Amélia Maio. 
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atravessado pela luz tal que o número de fotões 
produzido inicialmente (n,) é reduzido de 


o 


c 


— (n(L.) = ). 

Esse estudo é feito para fibras de diferentes diâme- 
tros e para diferentes dopantes. E do nosso conhecimento 
que as fibras se degradam sob o efeito de elevadas do- 
ses de radiação nuclear, e dado que nas situações físicas 
em que vão ser utilizados essas doses serão elevadas ( - 25 
krad/h), é do maior interesse o estudo detalhado da re- 
sistência destas fibras à radiação nuclear (y, n, e”) me- 
dindo o rendimento fotónico e o comprimento de ate- 
nuação antes e depois da irradicação das fibras. 

Por razões instrumentais é, em geral, necessário in- 
troduzir um guia de luz (material inactivo do ponto de 
vista de cintilação-plexiglas, por exemplo) entre a fibra 
cintilante e o fotodetector. Algumas técnicas foram en- 
saiadas para resolver este problema. Em colaboração com 
a Universidade de Neuchâtel (Suíça) e o CERN, estuda-se 
um novo método que consiste em irradiar parte da fi- 
bra activa com radiação U.V., destruindo a cintilação 
nessa zona. Deste modo a zona irradiada poderá fun- 
cionar como guia da luz produzida na zona da fibra não 
irradiada. 


Aparelhagem experimental e resultados 


A aparelhagem utilizada é típica em experiências de 
Física Nuclear. 

Como fonte de electrões (5 ) utiliza-se uma fonte ra- 
dioactiva de Sr”. A luz produzida pela passagem dos 
electrões que atravessam a fibra óptica de plástica é co- 
lectada na janela de um fotomultiplicador (PM2) de 
baixo ruído e boa resolução. A distância da fonte ra- 
dioactiva ao PM2 é variável permitindo estudar a trans- 
missão da luz produzida através da fibra. A saída do 
PM2 está colocada entre duas outras fibras Fl e F3 de 
igual diâmetro e que juntamente com a electrónica as- 
sociada servem de telescópio. Só são registados impul- 
sos no Multicanal quando há “simultaneamente”” (numa 
janela de 20 ns) detecção de luz nas fibras | e 3, i.e.,0 
sinal de coincidência entre a via 1 e a via 3 é utilizado 
como sinal de disparo na porta do multicanal. 

As propriedades das fibras são caracterizadas pelo 
número médio de fotoelectrões, fotões produzidos na fi- 
bra, ao ser atravessada por um electrão, transmitidos à 
distância /, e convertidos em electrões no fotocátodo. O 
número médio de fotoelectrões n,. é medido em função 
da distância / e o seu valor obtém-se a partir da análise 
do espectro registado no multicanal no MCA é propor- 
cional ao número de fotoelectrões produzidos em cada 
impulso). 

O espectro de fotoelectrões que se observa no multi- 
canal apresenta a forma gaussiana quando o número mé- 
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dio de fotoelectrões n,. é grande. No entanto, para fi- 
bras de diâmetro 4 < | mm, esse número é inferior a 10 
e daí os espectros apresentarem máximos e mínimos de 
distribuição mais complexa. De facto os espectros ob- 
servados ajustam-se à convolução de uma distribuição 
de Poisson cuja média é n, com uma distribuição de 
Gauss cuja largura é, essencialmente, função da estatis- 
tica de multiplicação do 1º dinodo do PM2. Na fig. 3 
apresenta-se um espectro-tipo para uma produção de um 
número médio de fotoelectrões n,.= 4. 

Como referimos anteriormente, foram estudadas fi- 
bras de diferentes características. Na fig. IV.4 estão re- 
presentados os resultados obtidos com fibras azuis e ver- 
des de diâmetro & = Imm e 0.5mm. | 


Estes resultados permitem-nos deduzir que: 


|) a produção de luz nas fibras de 0.5 mm compa- 
rada à das fibras de | mm é razoável sem que haja 
uma degradação no comprimento de atenuação; 
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Fig. IV.2 Dispositivo experimental para medição do nú- 
mero médio de fotoelectrões n,. e do comprimento de 
atenuação da fibra óptica cintilante F2. 
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Fig. IV.3 Espectro de fotoelectrões observados no se- 
lector de amplitude multicanal correspondente a um nú- 
mero médio de fotoelectrões n,> 4. 
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1) as fibras verdes, apresentam um comprimento 
de atenuação análogo ao das fibras azuis 
(Lo 15UM), 


No que respeita à resistência à radiação nuclear, re- 
sultados preliminares permitiram observar que, quando 
sujeitas a elevadas doses de radiação y (1 a 30 Mrad em 
3 a 10 dias) a degradação é mais pronunciada para as 
fibras azuis (Popop) do que para as fibras verdes (K27). 


Em) 


1,8 Eri) a 


Fig. IV.4 Número médio de fotoelectrões por electrão, 
ns em função da distância | a que a luz é transmitida. 


fibra azul de diâmetro 2? = Imm 
fibra verde de diâmetro & = Imm 

2 fibras azuis de diâmetro & = 0,5mm 
1 fibra azul de diâmetro q& = 0,5 mm 
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Fig. IV.5 Número de fotoelectrões por electrão atraves- 
sando uma fibra de Imm de diâmetro em função do 
tempo (t) de exposição a uma fonte intensa de UV. 


No desenvolvimento da técnica de utilização da pró- 
pria fibra cintilante como seu guia de luz, foram irra- 
diadas diversas fibras com radiação UV. e foi medido 
o número médio de fotoelectrões produzidos em função 
do número de horas de irradiação (t). A dose de irra- 
diação recebida pela fibra é proporcional a t. Na fig. 5 
apresentam-se os resultados obtidos com 6 lâmpadas Phi- 
lips TL de 40 Watt. No que respeita à degradação da 
produção de luz em função do tempo de exposição, re- 
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sultados análogos foram obtidos com fontes mais inten- 
sas (taxa de dose mais elevada) e menos intensas. 

Os resultados preliminares obtidos, permitem-nos 
prosseguir com o estudo da degradação das proprieda- 
des ópticas e de cintilação das fibras ópticas cintilantes 
de plástico, quando sujeitas à radiação nuclear ou à ra- 
diação U.V. embora por razões diversas como indicámos 
anteriormente. 
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RESUMO 


Analisam-se neste trabalho os processos básicos que 
ocorrem numa coluna de plasma criada e mantida por 
uma onda de superfície. Desenvolve-se um modelo teó- 
rico que permite determinar de forma consistente as dis- 
tribuições radiais da densidade do plasma e do campo 
eléctrico da onda, bem como a distribuição de densidade 
ao longo da coluna, em função das condições de funcio- 
namento. Mostra-se que existem leis de semelhança para 
este tipo de descarga que permitem a extensão dos resul- 
tados obtidos a condições experimentais diversificadas. 
O bom acordo obtido entre teoria e experiência permite 
concluir que o estudo destes plasmas abre o caminho para 
uma compreensão mais profunda das descargas HF. 


ABSTRACT 


This paper discusses the basic processes occurring in 
a plasma column sustained by a surface wave. The theo- 
retica! model here developed treats in a self-consistent way 
the radial distributions of the plasma density and of the 
SW field intensity and provides detailed predictions con- 
cerning the radial and the axial variations in the plasma 
properties as a function of the discharge parameters. H 
also predicts the existence of simple similarity laws which 
enable the extension of the results to different experimental 
circumstances. The theory is shown to agree well with the 
experiments. Flence, the study of these discharges opens 
the way to a deeper understanding of HF plasmas in ge- 
neral. 


|. INTRODUÇÃO 


Designam-se genericamente por ondas de superficie 
(O. S.) as ondas que se propagam ao longo da superfi- 
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cie de separação entre dois meios sem que haja conver- 
são de energia do campo associado com a onda para 
qualquer outro modo (sem radiação). A designação é 
apropriada já que, por um lado, a energia da onda é es- 
sencialmente transportada numa zona estreita junto da 
interface e que, por outro, a onda segue ao longo da in- 
terface, mesmo quando esta é encurvada, como se esti- 
vesse ligada a ela. 

Lord Rayleigh foi o primeiro a estabelecer que se po- 
dem excitar, através duma fonte de dimensões finitas, dois 
tipos de ondas num meio elástico semi-infinito: ondas 
espaciais (também designadas por ondas radiadas) que 
se propagam em todas as direcções; ondas de superfície 
que se propagam ao longo da fronteira. 

A situação é análoga no caso das ondas electromag- 
néticas. Se considerarmos, por exemplo, um cilindro fino, 
infinitamente longo, de material dieléctico sem perdas, 
imerso no ar, é possível a propagação duma onda elec- 
tromagnética de superfície ao longo da interface. A in- 
tensidade quer do campo eléctrico, quer do campo mag- 
nético, da onda decai radialmente, para uma distância 
r suficientemente grande, segundo a forma r-!? exp 
(-k, r), sendo k, a constante de decaimento. Na ausên- 
cia de perdas k, é real, pelo que não há variação de 
fase ao longo de r; não existe, portanto, fluxo de ener- 
gia perpendicularmente à interface, mas apenas ao longo 
dela. 

O decaimento do tipo exponencial da intensidade dos 
campos a partir da interface, para um ou ambos os la- 
dos, é uma característica das ondas de superfície. No caso 
dos dieléctricos usuais (isto é, excluindo os plasmas) só 
existe decaimento exponencial no dieléctrico de menor 
permitividade; no outro, a intensidade dos campos os- 
cila espacialmente (é, nomeadamente, o que acontece no 
interior do cilindro de dieléctrico no exemplo conside- 
rado acima). 

Contudo, se um dos meios for, por exemplo, um 
plasma frio e sem perdas, a sua permitividade 


e=1-05/w 


| 
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poderá ser negativa, o que se verifica quando a frequên- 
cia w for inferior à frequência plasma electrónica w,.. 
Trata-se do único dieléctrico conhecido com esta pro- 
priedade. Na superfície de separação de tal meio com 
um dieléctrico convencional de permitividade e >0 pode 
propagar-se uma O. S. com campos evanescentes em am- 
bos os lados da interface, desde que se verifique 
| €p | >e|9. O carácter evanescente puro apenas se ob- 
servará, no entanto, se ambos os meios forem semi-infi- 
nitos. No caso dum plasma de dimensões transversais 
finitas (por exemplo, coluna cilíndrica, fatia de plasma 
de geometria plana, etc.) o campo no plasma não é pro- 
priamente evanescente, embora a sua intensidade decresça 
da fronteira para o eixo. As O. S. com estas característi- 
cas chamam-se ondas de superfície puras*(> ? e serão 
as únicas consideradas neste artigo. 

O estudo das O. S. em plasmas, desde a sua desco- 
berta no trabalho pioneiro de Trivelpiece'? em 1958, 
tem motivado numerosos trabalhos teóricos e experimen- 
tais (veja-se o excelente artigo de revisão publicado re- 
centemente por Moisan et al)). As razões dum tal in- 
teresse foram, inicialmente, problemas de natureza 
fundamental em que intervêm as O. S., tais como o aque- 
cimento de plasmas por campos electromagnéticos, a ab- 
sorção de radiação laser por plasmas densos, as interac- 
ções não lineares como, por exemplo, as instabilidades 
paramétricas. Outro motivo de interesse residia no facto 
das O. S. poderem ser usadas para diagnóstico de plas- 
mas de baixa temperatura (por exemplo, plasmas de des- 
carga) permitindo medir a densidade, a frequência de co- 
lisão e a velocidade de deriva dos electrões t?). Em 
todos estes casos, no entanto, o plasma é suposto ser 
criado por mecanismos independentes da O. S., ou seja, 
esta propaga-se num plasma que se supõe pré-existente. 

Em anos mais recentes, no entanto, tem vindo a 
desenvolver-se um novo campo de investigação baseado 
na perspectiva inovadora da criação e manutenção de 
plasmas por ondas de superfíciel0%, Este tipo de plas- 
mas, que designaremos abreviadamente por plasmas de 
O. S., podem ser facilmente produzidos em laboratório, 
desde pressões da ordem dos militorr até algumas atmos- 
feras, usando quaisquer frequências desde alguns MHz 
até cerca de 10 GHz, o que permite cobrir todo o espec- 
tro desde as radiofrequências (rf) às micro-ondas. 

A densidade do plasma de O. S. pode variar tipica- 
mente entre 10!ºem-? e 10!cm-?, consoante a potência 
HF absorvida, a pressão do gás e o diâmetro do tubo 
de descarga. Por exemplo, pode-se produzir uma coluna 
de plasma com cerca de 2 m de comprimento, 2,5 cm 
de diâmetro e uma densidade média de cerca de 
10''cm-3, a uma pressão de alguns dezenas de militorr, 
usando 80 W de potência HF à frequência de 500 MHz. 
Os plasmas de O. S. são estáveis, calmos e perfeitamente 
reproductíveis, o que os torna substitutos interessantes 
dos plasmas de descarga convencionais (coluna positiva 
da descarga dc, descarga rf do tipo capacitivo, plasma 
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produzido numa cavidade micro-onda, etc.) para apli- 
cações de grande actualidade como o tratamento de su- 
perfícies, a gravura por plasma, a síntese química, a ex- 
citação de lasers, lâmpadas espectrais, etc. 

Para além do interesse destas aplicações, potenciais 
ou demonstradas, acresce que o estudo dos plasmas de 
O. S. tem contribuído significativamente para uma com- 
preensão mais aprofundada dos mecanismos básicos de 
funcionamento das descargas HF em geral (veja-se, por 
exemplo, sobre este assunto o artigo de Ferreira e Moi- 
sant!4). Tal facto deve-se, por um lado, à extrema ver- 
satilidade já mencionada de se poder varrer todo o es- 
pectro, desde as radiofrequências às micro-ondas, para 
se criar o plasma e, por outro, à possibilidade de rela- 
cionar localmente as propriedades do plasma com a po- 
tência absorvida (possibilidade que resulta do facto do 
plasma ser criado, neste caso, por uma onda progres- 
siva). Graças a estas características, a teoria dos plas- 
mas de O. S. atingiu já um grau de maturidade satisfa- 
tório, permitindo interpretar e explicar com boa precisão 
os resultados experimentais. 

Constitui objectivo central deste artigo discutir as pro- 
priedades dos plasmas de O. S. e os processos básicos 
que asseguram a sua manutenção. Por não terem ocor- 
rido, entretanto, progressos relevantes, parte do material 
científico aqui apresentado é, na sua essência, embora 
com variações de forma, o que se tomou como base para 
um curso dado sobre este assunto pelo autor no NATO 
Advanced Study Institute on Radiative Processes in Dis- 
charge Plasmas que teve lugar no Verão de 19859, 


2. PROPRIEDADES DAS ONDAS DE SUPERFÍCIE 
NUM PLASMA 


2.1] Equações dos Campos 


Considere-se uma O. S. propagando-se ao longo da 
interface (eixo -z) dum plasma com um dieléçtrico de 
permitividade e, e admitimos que os campos E e H va- 
riam segundo z e no tempo segundo a forma exp J (wt 
- KZ), em que x=8-ja, sendo 8=27/h o número de 
onda e q o coeficiente de atenuação. 

Descrevendo o plasma como um dialéctrico de per- 
mitividade «,, as amplitudes complexas dos campos 
obedecem às equações 


ada ia 


VxE=-jwn, H (1) 
V xH=jwee, E (2) 
V.H=0 (3) 
V.(e E)=0 (4) 


* Quando um dos meios em presença é um plasma com tp<O e existe um 
campo magnético aplicado paralelo à fronteira podem propagar-se igualmente 
ondas de superfície cuja intensidade do campo oscila espacialmente no inte- 
rior do plasma. 


2o 


A permitividade do plasma é, em geral, um número 


complexo. Está relacionada com a conductividade com- 
plexa do plasma, o, através da expressão 

= jd. (5 

E = ey AEE (5) 


Ú 


sendo o dado por!” 
4 So 42 
en. wu)? df 
E / ———— —— du 
m 5 v.+Jjw du 


Nestas expressões e e m representam a carga e a massa 
do electrão, e, é a constante dieléctrica do vácuo, n, é 
a densidade electrónica, v. é a frequência de colisão 
electrão-neutrão para a transferência de quantidade de 
movimento, u = mv?/2e é a energia do electrão expressa 
em volts e f(u) é a função de distribuição de energia 
(FDE) dos electrões com a normalização /. Prf(u) u!?2 
du = 1, Quando v. é independente de u a Eq. (6) reduz- 
-se à forma o = e2n./m (v. + jw); por outro lado, na au- 
sência de colisões (v.= 0) o é imaginário puro e a per- 
mitividade do plasma reduz-se à forma 


9=-— (6) 


(7) 


e,=1-(w/w), 


em que w.= (e“n,/em)”? é a frequência angular 
plasma electrónica. E » 
Tendo em atenção que = a -JK e, em que e, é 
o vector unitário do eixo Zz e Vo é a parte do opera- 
dor V que envolve unicamente as «derivadas em relação 
às coordenadas transversais, obtém-se facilmente das eqgs. 


(1) e (2) 


VixE-=-jwn, H (8) 
V.XE -juexE =-jomu, Ho (9) 
VixHr=jue e É, (10) 
VixH,-jxe,xH=jwe, e É, (11) 


em gue É, e H, designam as componentes transversais 
de E e H. E 

Multiplicando vectorialmente a eg. (9) por e, e com- 
binando com (11) obtém-se facilmente 


(12) 
E =-— o YE+ 2h (e xV,H), 
O cc É a 
c2 E, K aê En k 


em que c=(e, u,)!? é a velocidade da luz no vácuo. 
Procedendo de igual forma com as egs. (8) e (10) 
obtém-se agora 
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(13) 
— 1 K —+ 1 ty E. E —+ — 
Hj == sd Voy H,- E (e, x Vy E,) 
O é já MO ue 
a pk 7 K 


As egs. (12) e (13) constituem expressões gerais, re- 
lacionando as componentes transversais E, e H, com 
os gradientes transversos das componentes longitudinais, 
para qualquer onda progressiva. Reduzem-se a formas 
mais simples no caso de ondas TE (E, = 0), ou de on- 
das TM (H, = 0). Neste último caso, que nos interessará 
particularmente no seguimento, obtém-se simplesmente 


E,=- = Ve E, (14) 
it” 
2 

H = (e, x Er). (15) 


No caso dum plasma homogéneo (e, = const.) as 
equações que permitem obter E, e H, são particular- 
mente simples. 

Tendo em conta a identidade V x (V XA)= 
=" (V. AJ e A, obtém-se facilmente a partir das 


eqs. (1)-(3) 


v2 Erx(L e -x) E=0 (16) 
=» 2 oh | 
vi Er e) H=0 (17) 
portanto, E, e H, satisfazem às equações 
2 
vê E + e,-82) E,=0 (18) 
2 H Ea ad = > 
Ve H + e,-82) H,=0. (19) 


A hipótese de homogeneidade do plasma é, no en- 
tanto, inconsiste na generalidade dos casos e conduz a 
erros significativos. Por exemplo, se considerarmos uma 
coluna de plasma criada e mantida por uma O. S. em 
condições de atenuação fraca (w>>v,), é razoável des- 
prezar a inomogeneidade segundo o eixo do tubo, mas 
há que ter em conta, de forma consistente, a inomoge- 
neidade nas direcções transversais. Neste caso, partindo 
da Ea. (4) obtém-se VB=-E: Ve »/€p» pelo que a 
Eq. (16) adquire a forma mais geral 


+ 


+ 2 —a “+ —+ 
wz Er e) E+V (E,. 16,0 (20) 


p 
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Para uma onda TM E, é dado pela Eq. (14), a qual 
combinada com a Eq. (20) permite obter, neste caso, a 
seguinte expressão para E, 


De forma análoga, na Eq. (17) aparece um termo adi- 
cional escrevendo-se, agora, esta na forma 


” 


— 4 ] Quo "y 1 — TT 
V 5 dba RJ HA] o € Yy exE=0 


(22) 


Esta equação, conjuntamente com a Eq. (12), deter- 
mina H.. 


2.2 Ondas TM propagando-se ao longo dum plasma 
cilíndrico 


As ondas de superfície podem expressar-se, em ge- 
ral, como sobreposição de modos TE e TM. Contudo, 
é conhecido que apenas os modos TM podem satisfa- 
zer às condições de fronteira apropriadas, na interface 
plasma-dieléctrico, no caso de ondas com simetria axial. 
Assim, por exemplo, em geometria cilíndrica o único 
modo axialmente simétrico é o modo TM com m=0, 
em que m é o número de onda correspondente ao ân- 
gulo azimutal à; neste caso só existem as componentes 
E,, E, e H,. Todos os modos com m > são híbridos, 
possuindo simultaneamente componentes E e H, não 
nulas. 

As equações dos campos do modo TM no plasma 
obtêm-se facilmente aplicando as equações gerais (14), 
(15) e (21): 


SE, gifs K? A de dE, , 
dr? r w x € df dr 
Ei P 
o P 
(23) 


d E 
E, =-j) —— (24) 
w » dr 
3 EG-K 
w €, é 
H=2>.E, (25) 


Este modo é facilmente excitado usando um dispo- 
sitivo HF especial designado por “Surfatron”*(19:!), Tem 
sido este o modo usado, em geral, para a criação de plas- 


a a a a 


HF Power 


E 


Pyrex glass(eg) 
Plasma(eg) 


Fig. 1 Diagrama esquemático do dispositivo HF e do 
acoplador num Surfatron e da geometria do tubo de des- 
carga 


mas de O. S., pelo que nos limitaremos a este caso no 
seguimento. 

A Fig. 1 representa esquematicamente uma instala- 
ção experimental usada para produzir um plasma num 
tubo de vidro pyrex através duma O. S. “lançada” por 
um “'Surfatron”. 

As caracteristicas de dispersão das O. S. numa tal geo- 
metria obtêm-se a partir da continuidade das compo- 
nentes E, e H, nas interfaces plasma-vidro e vidro-ar 
(isto é, parar=a er=b, sendo a e b os raios interior 
e exterior do tubo, respectivamente). Para tal, basta ter 
em conta as Egs. (23)-(25) do campo no plasma e equa- 
ções similares quer no vidro, quer no ar exterior (com 
e, substituído por €,, constante, no vidro e por e, =, 
no ar). 

O cálculo das características de dispersão só é sim- 
ples, no entanto, no caso de ondas fracamente atenua- 
das («e<<), o que permite usar álgebra real, e admitindo 
que o plasma é homogêneo!!! !7. 1, Neste caso, as com- 
ponentes E, e E, no plasma são dadas pelas equações 
simples. 


E, (1)=E(o) 1, (Tr), (26) 


E (1) =E(0) 1, (W'r), (27) 
em que 1, e 1, são as funções de Bessel modificadas de 
ordem zero e um, respectivamente, E(o) é a intensidade 
do campo no eixo, e T é definido por 


- A 
2 


C 


(28) 
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ARGON, | = 360 MHz 
Tube: 2Za=25 mm 
:2b= 290.8 mm 
É q=4.52 
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Fig. 2 Características de dispersão duma onda de su- 
perfície com f = 360 MHz num plasma de argon 


A Fig. 2 mostra caracteristicas de dispersão calcula- 
das segundo este modelo simplificado e resultados ex- 
perimentais obtidos (experimentalmente estes resultados 
são obtidos. com w fixo e fazendo variar w,., pelo que 
o diagrama da figura se designa geralmente por diagrama 
de fase). Apesar das hipóteses simplificadoras usadas no 
cálculo, o acordo entre teoria e experiência é plenamente 
satisfatório, excepto para Sa elevado (o que é devido à 
não se ter considerado a inomogeneidade radial do 
plasma, como veremos em breve). 

Observam-se os seguintes limites no diagrama de fase. 
Quando w/w .—0 verifica-se que B-(w/c) Ve, em que 
|<e<e,. Significa isto que, quando w .>>w, o número 
de onda da O.5S. é próximo do que se observaria para 
uma onda plana electromagnética propagando-se num 
dieléctrico de permitividade e. No limite oposto, 58-00, 
verifica-se que w/w (1 + e) 12, ou seja que |€, | e, 
Este limite corresponde ao corte: a O. S. não se pode 
propagar se a densidade do plasma for inferior a um 
certo valor crítico, n,, dado pela relação |e | =e,. 
Numericamente, tem-se 

n. (cm?) = 1,245 x 10º fº(1 + e), (29) 
com f expresso em Hertz. 

A Fig. 3 mostra os valores de Ta e de Ba, calculados 
segundo este modelo, em função de w,. a/c, para vários 
valores da frequência e dos parâmetros geométricos a 
e b. Na vizinhança do corte tem-se T= 8 e Ta torna-se, 
por isso, elevado. 

As Eqgs. (26) e (27) revelam então que a intensidade 
do campo é muito maior junto da interface do que na 
zona axial do tubo. No limite oposto (w >>w, isto é, 
densidade do plasma elevada) verifica-se uma situação 
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Fig. 3 Curvas de fa (tracejado) e Ta (cheio) vs. 
We A/C num plasma de argon homogéneo: a = 0,45 cm, 
b=0,60cm, f=210MHz (A) e 2450MHz (B); 
a = 1,25cm, b = 1,49 cm, f = 210 MHz (D), 600 MHz (E), 
e 2450 MHz (F) 


semelhante, já que, nesse caso, se tem B<<T e Tau. 
a/c>>l. 

É importante notar que, quando a densidade do 
plasma é suficientemente elevada (0,.>>0), o quociente 
E (a)/E (o) é tanto maior quanto maior for a densidade 
do plasma; isto é, o campo concentra-se junto da inter- 
face plasma-vidro à medida que a densidade aumenta. 
Esta conclusão é igualmente válida para um plasma ra- 
dialmente inomogéneo (o que conta nesse caso é o va- 
lor da densidade média numa fatia de plasma) como ve- 
remos no seguimento. 


ELECTRIC FIELD INTENSITY 


Wu * UOlISOd jPIpEM 
ZH 009 


s=00/ "dm 


o: 
Õ 


Fig. 4 Distribuição radial das componentes radial 
(cheio) e axial (tracejado) do campo eléctrico duma O. S. 
TM propagando-se ao longo duma coluna cilíndrica de 
plasma homogéneo [da ref.] 
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A Fig. 4 mostra o perfil radial das componentes E, 
e E, no plasma, no vidro e no ar exterior, no caso 
f = 600 MHz e Woe/ 0 = 5. 

O efeito da inomogeneidade radial do plasma sobre 
o diagrama de fase da O. S. tem sido investigado por di- 
versos autores, usando a aproximação quase-estática 
(V xE=0>E=-VF). Trivelpiecel? e Ihé29 partiram 
da hipótese dum perfil parabólico da densidade, do tipo 


n(D=n, (O) [-K (DI, (30) 


em que O<K<1 constitui um parâmetro de forma do 
perfil. Os diagramas de fase w/w,. vs. Ba obtidas por 
ic, para diversos valores de K, são representados na 
Fig. 5 [(w. corresponde à densidade média n, = n, (0) 
(1-K/2)]. Verifica-se que para Ba<1 o efeito da inomo- 
geneidade é desprezável, tudo se passando como se a 
O. S. se propagasse num plasma homogéneo com den- 
sidade n.. Contudo, há diferenças significativas entre os 
dois casos quando Ba>1. A razão disso é, essencialmente, 
o facto do campo se concentrar, nesse caso, junto da in- 
terface, pelo que o corte fica determinado pela densi- 
dade na fronteira n. (a) e não pela densidade média (ou 
seja, w/w(a)+(1 + e)!* quando Ba-00). Sendo as- 
sim, obtém-se w/W al e MI-K)/Q-K/2)]* 
quando Sa-»0, valor que é sempre inferior ao que cor- 
responde a um plasma homogéneo com o mesmo valor 
de w 

Nos diagramas da Fig. 5 aparecem zonas de propa- 
gação retrógrada que podem, no entanto, como mostrou 
ic, ser eliminadas considerando ainda o efeito 
duma bainha junto da parede do tubo. Ilié apenas con- 


Fig. 5 feito dum perfil de densidade parabólico 
n,=n(0)[1-K (r/ a)*] sobre as características de disper- 
são do modo TM, para k = O(A); 0,4(B); 0,6(0); 0,8(D); 
e 0,9(E) [da ref.20] 
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siderou, contudo, o caso duma bainha sem electrões e 
dum plasma de densidade acima do valor crítico na fron- 
teira plasma-bainha, não tendo, portanto, analisado os 
efeitos da ressonância que ocorre quando w,. = w. Ora, 
se a densidade dos electrões decresce regularmente na 
bainha até se anular na parede, uma tal ressonância de- 
verá ocorrer algures, o que produzirá amortecimento dos 
modos de O. S. por excitação de ondas de Langmuir na 
região de ressonância'?), Este amortecimento virá adi- 
cionar-se ao produzido pela colisões electrão-neutro, co- 
lisões devido às quais os electrões do plasma extraem 
energia da onda. No momento actual, a teoria dos plas- 
mas de O. S. apenas tem em conta o amortecimento co- 
lisional, o que se tem revelado suficiente para a inter- 
pretação da generalidade dos resultados experimentais. 
A potência extraída da onda por uma unidade de 
comprimento, ao longo da direcção de propagação, de- 
vido às colisões elásticas electrão-neutro, é dada por 


E abs sd [ Roo) E? 27xr dr, 
PA Ros 


em que R,(o) é a parte real da conductividade, dada 
pela Eq. (6), e E? = E” + E?. Em regime fracamente co- 
lisional (atenuação fraca: v<<w) a Eq. (6) pode-se es- 
crever aproximadamente na forma 
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' en. . 
Re (0) = —— 
m 


a do (32) 


em que se definiu a frequência de colisão efectiva 


So 


v. = e 32 ve du. (33) 
Obtém-se assim 
dP a É 
Ei = e / v. (3 27mr dr, (34) 
o 


expressão que, no caso dum plasma homogéneo, se po- 
derá escrever ainda na forma 


d Pd =" Us 


dz 2 e) 


d 
(Lrey / Ei 2 dr. 
ni o 


O coeficiente de atenuação (colisional) da onda pode 
então calcular-se a partir da expressão (22) 


| dP./dz 
q=— abs” = a 
2 Va W 88) 
em que v, designa a velocidade de grupo e w, a média 
no tempo “da energia total (eléctrica + magnética + ciné- 
tica das partículas carregadas do plasma) por unidade 
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de comprimento, na direcção da propagação. Num 
plasma frio e em condições de fraca atenuação, Ve W, é 
igual ao fluxo total de energia electromagnética através 
dum plano perpendicular à direcção de propagação, ou 
seja, igual ao fluxo total da parte real do vector de Poyn- 
ting complexo E x H*/2. 


10 Ee 


creu) Telm! 


Fig. 6 Características de atenuação duma O. S. TM 
numa coluna cilíndrica de plasma de argon homogéneo: 
a=0,45cm, b= 0,60cm, f= 210 MHz (A) e 2450 MHz 
(B); a=0,30cm, b=0,45cm, f=2450MHz (C); 
a=1,25cm, b = 1,49 cm, f = 210 MHz (D), 600 MHz (E) 
e 2450 MHz (F) 


A Fig. 6 mostra valores de aw/v. em função de 
We! 4% calculados para um plasma homogéneo a partir 
do teorema de Poynting. Estes cálculos simplificados, que 
envolvem unicamente álgebra real, estão em excelente 
acordo com os resultados mais precisos que se podem 
obter a partir das soluções da equação de dispersão em 
álgebra complexa!?», É de notar, a partir dos resultados 
da Fig. 6, que a atenuação é pequena quando We >W 
(densidades do plasma elevadas), mas que ela se torna 
extremamente elevada na vizinhança da densidade de 
corte (convém ter em atenção que, na situação experi- 
mental que temos vindo a analisar, a densidade de corte 
corresponde aproximadamente ao fim da coluna de 
plasma; a densidade do plasma decresce regularmente 
ao longo da direcção de propagação, desde a saída do 
surfatron até se atingir a condição de corte, o que de- 
termina o fim da coluna). 


3. PROPRIEDADES DO PLASMA MANTIDO 
POR UMA ONDA DE SUPERFICIE 


3.1 Processos básicos 


Os processos básicos de manutenção dum plasma por 
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um campo eléctrico HF, em regime de absorção colisio- 
nal, são simples de entender. 

Devido às colisões electrão-neutro, o movimento dos 
electrões tem uma componente em fase com o campo, 
pelo que os electrões absorvem energia deste. Esta ener- 
gia é dissipada principalmente em colisões elásticas e ine- 
lásticas com os átomos do gás, incluindo-se nas colisões 
inelásticas a ionização dos átomos por impacto electro- 
nico com produção de novos electrões livres. Para que 
o plasma possa, então, ser mantido pelo campo HF é 
necessário que a taxa de produção de novos electrões por 
ionização compense exactamente a taxa de perda destas 
partículas por difusão para as paredes e, eventualmente, 
por recombinação. É necessário ainda, simultaneamente, 
que a potência absorvida por electrão, em média no vo- 
lume do plasma, seja igual à potência média perdida por 
electrão devido aos processos colisionais e de transporte 
para a parede circundante. 

A primeira destas condições pode ser adequadamente 
formulada através dum modelo do tipo hidrodinâmico, 
envolvendo as equações de continuidade e de transporte 
de quantidade de movimento, na presença do campo am- 
bipolar, para os electrões e os iões. A segunda pode ser 
formulada em termos cinéticos, isto é, através da pró- 
pria equação de Boltzmann, se admitirmos que os elec- 
trões se encontram localmente em equilibrio com o 
campo eléctrico HF e que os gradientes de difusão são 
suficientemente pequenos para não afectarem a distri- 
buição de energia dos electrões. Começaremos por ana- 
lisar a equação de Boltzmann e o balanço local de ener- 
gia. 


3.2 Função de distribuição de energia (FDE) 
dos electrões. Balanço local de energia 


A situação que analisaremos aqui é a dum gás de 
electrões na presença dum campo HF, E-=E, exp 
(ut)2%, Admitiremos, por hipótese, que: a) A ampli- 
tude do campo e os gradientes de difusão são suficien- 
temente pequenos para que a função de distribuição de 
velocidades dos electrões possa ser representada por um 
desenvolvimento em harmónicas esféricas limitado aos 
dois primeiros termos (parte isotrópica + uma anisotro- 
pia de 1? ordem). Tal hipótese requer, em particular, que 
a secção eficaz total de difusão seja muito maior que 
a de colisão inelástica e que o livre percurso médio dos 
electrões seja muito inferior a qualquer das dimensões 
do recipiente; b) A frequência w é suficientemente ele- 
vada para que a parte isotrópica da distribuição não se 
modifique de forma apreciável durante um ciclo de os- 
cilação do campo. Esta hipótese é válida desde que 
w>>ô,!, em que ô, é o tempo característico de relaxa- 
ção da energia por colisões com os átomos (note-se que, 
para colisões elásticas, ó! = (m/M)v., sendo M a 
massa dos átomos, enquanto que, para as colisões ine- 
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lásticas, 6,! = v,, sendo v, a frequência de colisão ine- 
lástica total): o Os fluxos no espaço de configuração e 
no de velocidades resultantes dos gradientes de difusão 
e do campo de carga de espaço são pequenos e podem 
ser desprezados, isto é, a distribuição de energia dos elec- 
trões é predominantemente determinada pelo campo HF 
aplicado e pelas colisões elásticas e inelásticas (uma dis- 
cussão dos efeitos do campo de carga de espaço sobre 
a distribuição pode encontrar-se em. 26), 

Com estas hipóteses e seguindo o desenvolvimento 
feito por Holstein??, a parte isotrópica da distribuição 
é dada pela equação!” (equação de Boltzmann homo- 
génea) 


df 


- de [3/* U. Qu, SE + f)] = v(u) f(u) Vu- 
37) 
-Lvu(u+V)f(u+V) V(u+Vo, 
ri 
em que 
at A K 
u(W=— dra 2 (38) 


representa a energia média extraída pelos electrões de 
energia u do campo aplicado, por colisão (portanto, o 
produto v, u. é a potência transferida). Na Eq. (37) v 
representa a frequência de excitação do nível j, com li- 
miar V,, a partir do fundamental e v,=L v;; a soma L 

J | 
inclui a ionização, desde que se despreze a energia do 
electrão secundário que nasceu; a normalização de f é, 
como anteriormente, de fvu du=1. 


Multiplicando ambos os membros da Eq. (37) pela 
energia do electrão e integrando sobre todas as energias, 
O<u< co, obtém-se a equação de conservação da energia 


2m 


E? 
(e n,)! R.(0) = AR cy u>+& Vi<up, (39) 


em que o é a conductividade, dada pela Eq. (6), e <> 
representa a média duma dada quantidade em toda a 
gama de energia. A Eq. (39) exprime simplesmente o 
facto que a energia absorvida do campo é igual à que 
é perdida por colisões elásticas (1º termo do 2º mem- 
bro) e inelásticas (2º termo), por unidade de tempo e 
em média por electrão. 

A Eg. (37) mostra que a FDE só depende das variá- 
veis exteriores E, w e N (densidade do gás) através da 
quantidade U.. Quando v.(u)>>w, para qualquer ener- 
gia u, ter-se-á u,= e E2/mv) a(E/N), em que 
E=E,/V2éo valor eficão do campo; a FDE será pre- 
cisamente a mesma que a que corresponderia a um 
campo d.c. de amplitude E. O único parâmetro indepen- 
dente é, neste caso, E/N. No limite oposto, w>>v. (u), 


ter-se-á u, - e E?/mwa(E/w) e a FDE só depende 
do parâmetro E/w Em qualquer situação intermédia a 
FDE depende simultaneamente de E/N e de E/w, ou, 
se preferirmos, de E/N (ou E/w) e de N/qw Como 
v.(u)a N, é preferível usar em vez de N/w, que tem di- 
mensões, um parâmetro adimensional do tipo v,../u, 
sendo v.. a frequência de colisão para os electrões 
duma dada energia (a qual poderá ser escolhida arbi- 
trariamente). 

Quando v. é independente da energia do electrão é 
possível reduzir a dependência da FDE a um único pa- 
râmetro, para qualquer valor de v./w, introduzindo um 
campo efectivo!?8. 29) 


.=E aaa To (40) 


Neste caso, u. = (e/m) (E,/v) a (E /N), pelo que 
a FDE só depende, de facto, de E,/N para qualquer va- 
lor de v./«, 

Por outro lado, o balanço da potência, expresso pela 
Eq. (39), pode escrever-se neste caso (em eV's'!). 


cE*  2m 
UV; Uç = | = —— 
mM v, 


<u Geo da Vi<u> (41) 

O conceito de campo efectivo é poderoso (todas as 
grandezas médias referentes aos electrões só dependeriam 
de E./N) mas apenas interessa para modelos físicos 
simplificados, já que, em geral e com excepção do hélio 
e do hidrogénio, v. depende fortemente da energia do 
electrão. 

No caso mais geral em que v, depende de u é inte- 
ressante classificar os diferentes regimes de funciona- 
mento através dum diagrama de E/N vs. E/«w, como se 
mostra na Fig. 7 para o caso do argon'?. A linha recta 


- 
E 
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Fig. 7 Diagrama de E/N vs. E/w: (A) R=w/vy= 1; 
(B) R'=w/v,=10; (€C) E/w=3x 10º Vem''s. A 
curva a traço-ponto é o limite assimptótico da linha 
R' = 10 para pequenos valores de E/w 
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A corresponde à condição R=w/vym=1, em que 
vm/N=3.2x 107º cm */s constitui o valor máximo de 
v./N (que ocorre para u = 12 eV). No semi-plano supe- 
rior (R>1) a FDE torna-se unicamente função de E/w 
(regime HF), enquanto que no inferior (R<1) será uni- 
camente função de E/N, exceptuando obviamente uma 
zona de transição entre as duas regiões. A condição 
w>v.m poderá ser, todavia, demasiado severa para o es- 
tabelecimento do regime HF, e poder-se-ia substituir, por 
exemplo, por uma condição do tipo w>v., em que v, é 
definido pela Eq. (33). A curva B corresponde a 
w=10v., o regime HF estará certamente estabelecido 
no semi-plano acima desta curva. Finalmente, a linha 
vertical corresponde ao valor de E/w para o qual a perda 
de energia por colisão elástica e inelástica é precisamente 
igual, em regime HF; para valores de E/w inferiores 
(zona à esquerda) a perda elástica predomina e o corpo 
da distribuição torna-se Maxwelliano com temperatura 

=(eM/3m?) (E/w)(l6. 24); indicam-se no diagrama 
os valores de T, em correspondência com os de E/q, 
nesta situação. 
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Fig. 8 Função de distribuição de energia dos electrões 
no argon para E/N = 6,5x 10-1º Vem? e R=w/vy=O 
(A); 1,25 (B); 9,35 (C); 15,6 (D); 23,4 (F); 31,2 (G) 


A forma da FDE modifica-se substancialmente à me- 
dida que se passa da situação w<<v., para a de w>>v.. 
A Fig. 8 mostra a FDE no argon para um valor cons- 
tante E/N=6,5x 10-16 Vem” e vários valores de 
R =w/v. desde zero (DC) a 31,2 (HF). À medida que 
R aumenta o ponto de funcionamento no diagrama da 
Fig. 7 desloca-se horizontalmente da direita para a es- 
querda. desde a região DC para a região HF. O corpo 
da FDE torna-se uma função rapidamente decrescente 
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Fig. 9 Função de distribuição de energia dos electrões 
no argon para uma mesma energia média de 3,5 eV e 
os seguintes valores de w/v,. = O (A); 0,5 (B); 0,8 (C); 
>> 1 (D) 


da energia, o que se traduz por um aumento significa- 
tivo da fracção de electrões com energias elevadas (cauda 
da FDE). A razão física disso é fácil de entender: a po- 
tência transferida v. u. é máxima quando v.(u)=w, 
logo, sendo v.(u) uma função crescente até u = 12 eV, a 
transferência de potência em regime HF é maior para 
electrões com energias da ordem da dezena de electrões- 
-volt, o que enriquece substancialmente a cauda da dis- 
tribuição. 

Este enriquecimento relativo da cauda é bem demons- 
trado na Fig. 9 em que se mostram várias distribuições 
com a mesma energia média, nos casos DC, HF, e in- 
termédio. Uma consequência importante deste facto é 
que os processos de excitação de limiar mais elevado e 
a ionização se tornam mais eficazes em regime HF, como 
se ilustra na Fig. 10. Nesta figura mostra-se a fracção 
de potência perdida pelos diversos canais colisionais nos 
casos DC e HF, em função da potência total dissipada 
por colisões, em média por electrão, para uma densidade 
unitária do gás, 0/N*; as curvas correspondentes ao caso 
w>>v, estão transladadas para menores valores de O/N 
em relação às do caso w<<v.. 

Como se pode verificar, para um dado 9/N a ftrac- 
ção de potência gasta na ionização é maior no caso HF 
que no DC; expressando isto noutra forma, para man- 
ter uma descarga compensando as perdas de electrões 


* Note-se que 


2m U.. u 
<ul> + EV<-ls, 
M N j 


= (en)! Re() E? = 
e) q? 
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Fractional power transfer(%) 
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Fig. 10 Perdas de energia relativa dos electrões no ar- 
gon pelos diferentes canais colisionais em funções da po- 
tência média absorvida por electrão para densidade uni- 
tária do gás, nos limites w>>v,, (a cheio) e w<<v,, (a 
tracejado): E — colisões elásticas, M — 3p, +ºp,; R — 
'p; P — %p,; F — transições proibidas; H — estados 
permitidos elevados; I — ionização 


pela ionização em volume, é necessária uma menor po- 
tência média absorvida por electrão no caso HF do que 
no caso DCO9. Para uma mesma potência absorvida 
por unidade de volume o plasma HF terá, portanto, 
maior densidade que o plasma DC, assim como uma 
maior população de estados de limiar de excitação ele- 
vado. Estas conclusões concordam bem com diversos fac- 
tos experimentais referidos na literaturaf». 


Para o balanço de criação-perda dos pares electrão- 
-jão (problema que, como referimos, é fundamental no 
estudo do equilíbrio do plasma) é necessário conhecer 
o coeficiente de ionização por colisão electrónica 
C,=<o, v> em que o, é a secção eficaz de ionização. 
Em princípio C,, já que envolve uma média sobre a 
FDE, deverá exibir um tipo de dependência em relação 
a E/N e E/w (ou E/N e v../w) como o da própria 
FDE. No argon, contudo, Ferreira e Loureiro? mostra- 
ram que é possível definir um campo efectivo para a io- 
nização da forma (note-se a analogia com o campo efec- 
tivo definido no caso v.(u) = const.) 


(42) 


é (v2 + w2)1/2 


com v.,./N=2x 107 cm*/s, o qual permite expressar C, 
como função da única variável Ei/N, para qualquer 
frequência (isto é, qualquer valor de w/v..). A Fig. 11 
mostra o coeficiente de ionização no argon, obtido a par- 
tir da resolução da equação de Boltzmannº), 


3.3 Balanço global de energia 


O balanço global de energia duma descarga HF pode 
escrever-se na format!s 32) 


P, A(n)=0 V,n, 


em que P, é a potência HF incidente no plasma, A (n,) 
é a fracção dessa potência que é absorvida no volume 
do plasma V,, n. é a densidade média em V,, e O re- 
presenta a potência média perdida por electrão em coli- 
sões de todas as espécies no volume V.. Em princípio 
O depende explicitamente de quantidades como a pres- 
são do gás, as dimensões do recipiente, a densidade mé- 
dia n., a frequência de funcionamento e ainda, em 
princípio, das distribuições espaciais da densidade do 
plasma e do campo HF no interior do plasma. Esta 
quantidade desempenha um papel fundamental na teo- 
ria dos plasmas HF pelo que, no seguimento, a analisa- 
remos em pormenor no caso particular dos plasmas de 
O. 8. 

O termo no primeiro membro da Eq. (42) representa, 
por assim dizer, o comportamento electrodinâmico do 
aplicador quando tem por carga o próprio plasma; deste 
ponto de vista, o plasma é fenomenologicamente repre- 
sentado como um dieléctrico cuja permitividade, para 
uma dada frequência, é unicamente determinada pela sua 
densidade média e pela frequência de colisão. A frac- 
ção de potência incidente que é absorvida, A (n,), de- 
pende, naturalmente, do tipo de aplicador e de instala- 
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Fig. 11 Coeficiente de ionização no argon em função 
do campo efectivo reduzido para vários valores de w/ Voe 
entre zero e infinito (pontos). A curva representa a me- 
lhor interpolação destes pontos 
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ção. Assim, por exemplo, no caso de plasmas produzi- 
dos dentro duma cavidade ressonante, A (n,) tem a 
forma duma curva de ressonância cuja amplitude de- 
pende do acoplamento entre o excitador e a cavidade. 
No caso dos plasmas de O. S., ou de onda progres- 
siva em geral, A (n,) pode ser resolvido espacialmente. 
Com efeito, designando por P (z) o fluxo total de po- 
tência através do plano z = const., a potência absorvida 
numa fatia de plasma entre z e z + dz pode escrever-se 
P(z) dA(n)=2o(n,) P(z) dz, (43) 
em que q é o coeficiente de atenuação, anteriormente 
definido pela Eq. (36). Tendo em conta a Eg. (42), o ba- 
lanço global de potência na fatia de plasma poderá, pois, 
escrever-se na forma 
2a (n.) P(z)= On, 7 a?. (44) 
Mas, a potência absorvida por unidade de volume pe- 
los electrões é dada por dP,. /dV = R,(o) E”, o que se 
pode ainda escrever, no caso de v. ser independente de 
u, 


AP a E É (45) 
dV mo uv. 


C 


sendo E, o campo efectivo”, Integrando a Eq. (45) na 
secção transversal e combinando com (44) obtém-se fa- 
cilmente 


(46) 


Convém escrever esta equação duma forma um pouco 
diferente que será útil para pôr em relevo a existência 
de relações de semelhança. Assim, escrevê-la-emos na 
forma 


o E oO (Fey e. E* 2x dx, (47) 
N m (v./N) A n. 
em que E,, é o campo efectivo no eixo e E= E /E 


eo” 


Verifica-se, assim, que 0/N é uma função explícita 
de E.,/N, o qual, como veremos, constitui um valor 
próprio do problema de valores na fronteira que resulta 
do equilibrio criação-perda dos pares electrão-ião. Por 
outro lado, O/N depende, através do termo integral, das 
distribuições espaciais da densidade e da intensidade do 
campo no plasma. Estas distribuições espaciais, tal como 
o valor de E,,/N, podem ser consistentemente determi- 
nados a partir da formulação do equilíbrio criação-perda 
acima referido. 

Quando v. (u) + const. os raciocínios são idênticos 
mas as expressões são mais complicadas, por não se po- 


der fazer pleno uso do conceito de campo efectivo. Só 
não será assim no caso v (u)<<w, ou seja, para fraca 
atenuação, vindo neste caso [ver Eqgs. (32)-(34)] 


E sql / Dep 
N m uu o n 


em que E, designa agora o valor eficaz da intensidade 
do campo no eixo e E=E/E.. 


Ccox dx, (48) 
N 


3.4 Equilíbrio de ionização-perda. 
Campo de manutenção 


Num plasma quase-neutro (comprimento de Debye, 
Ap, muito menor que qualquer dimensão do recipiente) 
o transporte dos electrões e dos 10es para a parede pode 
ser descrito pelo modelo hidrodinâmico de dois momen- 
tos, isto é, usando as equações de continuidade e de 
transporte de quantidade de movimento para os electrões 
e os i0es8%, Estas podem escrever-se 


V.(nv)=vn (49) 
cj Ei TV (Do (50) 
n 
(v. Dire Ó o Sl ven (51) 
M M ím 


em que n=n,=n, é a densidade do plasma, v é a ve- 
locidade média de aii comum às duas espécies 
(transporte ambipolar), E, é o campo de carga do es- 
paço, u. e D. são a mobilidade e o coeficiente de difu- 
são dos electrões, v =NC é a frequência de ionização, 
T, é a temperatura dos i iões, ev, é a frequência de co- 
lisão para a transferência de quantidade de movimento 
ião-neutro incluindo os efeitos da transferência de carga 
simétrica e das colisões elásticas. 

Este sistema de equações permite fazer surgir natu- 
ralmente a fronteira plasma-bainha de carga de espaço 
(junto da parede)º*”, aparecendo esta fronteira como os 
pontos em que a velocidade ambipolar v se torna igual 
à velocidade isotérmica iónica do som [K(T. + T)/M]'” 
(o que satisfaz à condição de Bohm). 

Note-se que T. representa, aqui, a temperatura dos 
electrões quando a FDE é Maxwelliana; quando assim 
não acontece T, deve ser substituída pela energia carac- 
terística u, = D./u 89. 

Como a bainha de carga de espaço tem espessura des- 
prezável (A,<<A, sendo A o comprimento de difusão 
característico do recipiente) pode-se admitir que a con- 
dição v=[K(T.+ T)/M]'"* ocorre na parede; desta 
forma, para uma dada geometria do recipiente, as Egs. 
(49)-(51 constituem um problema de valores na fronteira. 
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No caso dum plasma isotérmico (T. espacialmente 
constante), por exemplo, fixando a densidade do gás N 
e supondo T, conhecido, só existe solução para um va- 
lor de T, bem determinado (o qual determina o valor 
do coeficiente de ionização). A solução do problema 
pode então expressar-se sob a forma usual duma carac- 
terística de descarga de T, em função de NA. Contudo, 
uma tal característica carece de sentido quando os elec- 
trões não são Maxwellianos: neste caso, como o coefi- 
ciente de ionização C, e a energia característica u, são 
funções do campo reduzido (E/N no caso dc. ou 
w<<v.; E/w no caso w>>v.; E./N no caso de v (u) = 
= const., para qualquer w), as Egs. (49)-(51) determinam 
o campo reduzido em função de NA. 

Verifica-se, assim, que o campo eléctrico no interior 
do plasma não tem um valor arbitrário, mas sim, um 
valor bem determinado (chamado campo de manuten- 
ção), quando se fixam as condições de funcionamento 
(N, w, geometria, etc.). 

Esse valor é precisamente aquele para o qual a fre- 
quência de ionização compensa exactamente a taxa de 
perda dos electrões e dos 10es por transporte para a pa- 
rede, onde se vão recombinar. 

Se o plasma não for isotérmico, o que acontece no- 
meadamente no caso dos plasmas de O. S. visto que o 
campo não é homogéneo no interior do plasma, o pro- 
blema é um pouco mais complicado porque v, e u, va- 
riam espacialmente. Essa variação é determinada pela 
própria variação do campo, a qual, por sua vez, depende 
da distribuição espacial da densidade do plasma [ver, por 
exemplo, as Eqgs. (23) e (24)]. Para tratar o problema 
duma forma consistente é, por isso, necessário resolver 
simultaneamente as Eqgs. (23), (24) e (44)-(51), após ter 


p=0.1 torr; (= 600 MHz 


T27.5 x 10 em? 
--—— Tg = 10" cm? 


/ 
eme RE 2.5x10"enr3 f 


r/a 


Fig. 12 Distribuição radial de densidade electrónica 
num plasma de argon mantido por uma 0. S. TM, para 
p=0.1 Torr, f = 600 MHz e vários valores da densidade 
média 
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p=Ottorr; 1=600 MHz 


e 

e n 

a = 4x10 cm? 
c 


fes 2.5x10" em? 


N 


nes <10"em-* 


0.50 


0.25 


1.0 
r/a 
Fig. 13 Distribuições radiais das intensidades das com- 
ponentes axial e radial do campo eléctrico duma O. 8. 
TM propagando-se num plasma de argon inomogéneo, 
para diversos valores da densidade electrónica média 


0.5 


resolvido a própria equação de Boltzmann que analisá- 
mos no parágrafo 3.2. Note-se, no entanto, que a natu- 
reza do problema de valores na fronteira se mantém; 
neste caso, é agora o valor do campo reduzido no eixo 
do plasma que constitui um valor próprio do problema. 
Este campo de manutenção no eixo depende agora não 
só de N, w, À, como no caso isotérmico, mas também 
da própria densidade média do plasma n, (note-se que 
as Egs. (23) e (24) mostram que a distribuição espacial 
do campo depende não só da distribuição espacial rela- 
tiva da densidade mas também do seu valor absoluto, 
o qual pode ser medido, por exemplo, pelo seu valor no 
eixo ou pela média n,). 

Nas figuras 12 a 15 mostram-se resultados obtidos 
anteriormente pelo autorC”, através desta formulação 
consistente, para o caso dum plasma de argon mantido 
num tubo de pyrex (e, = 4.52) cilíndrico por uma O. 5. 
TM fracamente atenuada. A medida que aumenta a den- 
sidade média do plasma, o campo total (E?+ E?!” 
concentra-se cada vez mais junto da parede, o que au- 
menta consideravelmente os valores das frequências de 
ionização e de excitação por colisão electrónica nesta re- 
gião em relação aos valores no eixo. Uma maior taxa 
de ionização na periferia aumenta os gradientes de di- 
fusão nesta região ao mesmo tempo que os diminui na 
zona central, o que é bem evidenciado pelos perfis de 
densidade obtidos. 

Como se vê na Fig. 14 a inomogeneidade radial do 
plasma afecta consideravelmente a distribuição espacial 
do campo em relação à que se observaria se o plasma 
fosse homogéneo. Note-se ainda que a componente radial 
do campo é desprezável face à axial para n, bastante 
acima da densidade do corte (n,. = 2,47 x 10'º cm, 
para f=600 MHz e e, = 4,52). 
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f = 800 MH: 


p=0.1 torr ; 


E/E, 


— SW PLASMA 
-——— HOMOGENEOUS 


73 | PLASMA 


1.50 


1.25 


0.5 1.0 
ra 
Fig. 14 Distribuições radiais da intensidade do campo 
eléctrico duma O. S. TM num plasma inomogéneo (a 
cheio) e num plasma homogêneo (a tracejado), para a 
mesma densidade electrónica média 


304/m ; 600 MHz 
A=-po0os torr 
B-p=0.1 torr 
C-ps 0.2 torr 
D=p=z 0,4 torr 


Fig. 15 Distribuições radiais da densidade de átomos 
de argon excitados na configuração 3 p*6 d num plasma 
mantido por uma O. S. TM, para uma mesma potência 
absorvida por unidade de comprimento e para várias 
pressões 


A Fig. 15 mostra perfis radiais calculados de átomos 
excitados na configuração 3 p” 6d para várias condições 
de funcionamento. A forma destes perfis pode ser facil- 
mente compreendida atendendo a que a população des- 
tes estados é, localmente, proporcional ao produto n, 
C,, sendo C, o coeficiente de excitação electrónica cor- 
respondente. A densidade n. decresce com o raio, mas 
C, aumenta e tanto mais rapidamente quanto maior for 
n.; O produto n, €, poderá pois passar por um máximo 
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para certa posição radial, o qual será tanto mais acen- 
tuado quanto maior for n,. E precisamente o que ilus- 
tra a Fig. 15, já que o aumento de pressão, mantendo 
a potência absorvida por unidade de comprimento, se 
traduz por um aumento de n,. Estes perfis calculados 
estão inteiramente de acordo com observações experimen- 
tais usando técnicas de espectroscopia de emissão! *7, 


3.5 Leis de Semelhança 


O plasma de O. S. fracamente atenuada (w>>v,) 
obedece a leis de semelhança simples que permitem a 
extensão de resultados obtidos numa dada situação par- 
ticular a outras situações experimentais"”. Para prová- 
“lo, convém escrever as equações básicas numa forma 
apropriada. 


a) Equações de continuidade e de transporte de quan- 
tidade de movimento: 


1d ,y=-(Na)Cn (52) 
E OX 
DN e do 1 dn - 
Na KT. dx n dx 
prot SP Bh à 
dx M dx M mn dx 
(54) 


- (Na) (C, + E) v. 


b) Equações do campo eléctrico da O. S. no plasma: 


(55) 
dE AM B* 1 dn, dE 

— — + -— — * (Va) E,=0 

dx? Rise DP mn dx” dx ind 
Ba d E E 
E, =——— ————+&., 56 
"o (Ta) dx (2) 

c) Balanço local de energia dos electrões: 

(57) 

Dl ht TO 

N m Now M N A 


Nestas equações x=r/a é a coordenada radial nor- 
malizada, E, =E,/E, e E,=E,/E, são as intensidades 
das componentes do campo normalizadas em relação à 
intensidade no eixo, é (x) é o potencial do plasma, v, 
e T são as grandezas definidas nas Egs. (33) e (28), res- 
pectivamente. 

Como já vimos, quando >> verifica-se que 
"= w a/c>>B. Na prática este limite é logo atingido 
quando w,.. excede um pouco o valor de corte 


pe 
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Ure = (1 + e,)“2, pelo que esta aproximação é válida em 
quase toda a coluna de plasma (exceptuando apenas a 
região em que esta termina e, eventualmente, a região 
periférica junto da parede do tubo). 

Neste caso as Egs. (55) e (56) podem escrever-se apro- 
ximadamente 


2 | tw 
d E EM d E, (Cel y de É. = () (58) 
dx? X dx C n 
E<<E,. 28) 


Desta forma verifica-se que apenas existem dois pa- 
râmetros independentes no problema: Na e w a/c. Tal 
facto significa que as distribuições espaciais de densi- 
dade n.(x)/n., de intensidade do campo E (x) e de tem- 
peratura T. (x), assim como a intensidade do campo re- 
duzido no eixo E /w, são unicamente funções daqueles 
dois parâmetros. 

Sendo assim, imediatamente se conclui que a frequên- 
cia média de colisão para densidade de gás unitária, 
v./N, e a potência média absorvida por electrão numa 
fatia dz de plasma, para densidade de gás unitária, O/N 
[ver Eq. (48)], são unicamente funções de Na e de 
Wed/C: 

Cc -f (Na, Ce2) 
i c 


(60) 


N 
O we À 
ET fo; (Na, TR (61) 
Estas conclusões constituem leis de semelhança para 
este tipo de descarga. Verifica-se que estão plenamente 
de acordo com os resultados experimentais, como a se- 
guir se mostra para O caso do argon. 
A Fig. 16 mostra uma comparação dos valores me- 
dios radiais de v./p teóricos com dados obtidos de 
várias experiências”. 2) (p é a pressão para uma tempe- 


o 


e /p 10º torr=! s=1) 


Pe 


“Nydods Day 


100" 1 
pa itorr em] 


10-2 


Fig. 16 Frequência de colisão efectiva a pressão unitá- 
ria num plasma de argon mantido por uma O. S. TM 
em função do produto pa. Curva a cheio-teoria; pontos 
- resultados experimentais das refs.(”7. 2 


107!º 


0/p (Watt torr=") 
S 
1 


107! 1 
pa itorr cmi 


Fig. 17 Potência média absorvida por electrão a pres- 
são unitária requerida para manter uma coluna cilíndrica 
de plasma de argon em função do produto pa, nas se- 
guintes situações: i) Plasma isotérmico: w>>v,., FDE 
não Maxwelliana (A); w>>v.., FDE Maxwelliana (B); 
w<<v,. FDE não Maxwelliana (D); ii) Plasma mantido 
por O.S. TM, para w>>v,.: teoria com FDE Maxwel- 
liana (C); resultados experimentais de Dervisevic et alt” 
(pontos) 


ratura de 300k, isto é p=1Torr corresponde a 
N = 3,22 x 10!ºcm-?). O acordo entre teoria e experiên- 
cia é satisfatório, o que indica, em particular, que a onda 
é principalmente atenuada pelas colisões electrónicas à 
medida que se propaga. 

A Fig. 17 mostra, para fins de comparação, valores 
calculados de 9/p em função de pa para diversas situa- 
ções (casos DC e HF supondo o plasma isotérmico, 
plasma de O. S. com electrões Maxwellianos e plasma 
isotérmico com electrões Maxwellianos). No caso do 
plasma de O. S. seria de esperar que O/p dependesse 
também de w,. a/c, mas esta dependência é, de facto, 
tão ligeira que não se evidencia nos resultados da figura. 
Verifica-se que os cálculos de 0/pvs. pa no caso do 
plasma de O. S. com electrões Maxwellianos estão in- 
teiramente de acordo com os resultados experimentais 
obtidos por Dervisevic et alt”), autores que investigaram 
plasmas de O. S. mantidos em capilares, usando várias 
frequências na banda das micro-ondas. (Note-se que nes- 
tes plasmas n, é relativamente elevado, razão por que a 
FDE deve ser aproximadamente uma Maxwelliana). 

Finalmente, a Fig. 18 mostra o campo efectivo redu- 
zido no eixo E,,/N = (E,/N) (v,./w) em função de Na, 
para dois valores de w a/c admitindo que a FDE é 
Maxwelliana. Mostram-se, igualmente, para comparação, 
resultados obtidos para um plasma isotérmico nos ca- 
sos DC e HF com a FDE obtida da equação de Boltz- 
mann, assim como no caso HF com electrões Maxwel- 
lianos. 
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Ee/N(107'º v em?) 


9 “o 10 
Na(cm"?) 


Fig. 18 Campo efectivo reduzido no eixo necessário 
para manter um plasma cilíndrico de argon em função 
do produto Na para w>>v,. e nas seguintes situações: 
i) Plasma isotérmico: FDE Maxwelliana (A) e FDE não 
Maxwelliana (D); ii) Plasma mantido por O. S. TM, com 
FDE Maxvwelliana, para w, a/c=0,75 (B) e w 
a/c=1,5 (€) 


pe pe 


3.6 Propriedades axiais do plasma de onda de superfície 


Nos parágrafos anteriores analisaram-se as proprie- 
dades radiais do plasma de O. S. e o equilíbrio de ioni- 
zação-perda e de energia em cada fatia dz. É agora fácil 
descrever a evolução da densidade média n, ao longo 
do eixo em função das condições de funcionamento. 
Com efeito, a partir da Eq. (44) obtém-se 


20 (n)P(2) 


c O miz (62) 


logo, derivando em ordem a z esta equação e combi- 
nando de novo com (44), virá 


dn, Ba Edo n. da 


63 
dz a dn, 163) 


. 


Desprezou-se aqui a dependência de O em relação a n. 
por ser, como vimos, muito fraca. 

Usando as características de atenuação (Fig. 6) esta 
equação poderia ser facilmente integrada numericamente 
para obter n.(z) em cada caso. A densidade no início 
da coluna n,(0) (2=0 corresponde à saída do surfa- 
tron) é simplesmente dada pela solução da equação 


(64) 


em que P é a potência total incidente em z=0, O/N 
é determinado pelo valor de Na e «,=c [n.(0)). 

É, no entanto, útil evitar a integração numérica 
fazendo aproximações razoáveis que permitam obter 
expressões analíticas. Como se depreende da Fig. 6, log 
(ow/v.) varia quase linearmente com log (w,./w), ex- 
cepto na vizinhança do corte (termo da coluna), pelo que 
usaremos a aproximação 


- (65) 


gato) ps, 
E] 


E 


em que g (w, a) e h se obtêm facilmente, em cada caso, 
por interpolação dos resultados numéricos através da 
Eq. (65). A integração da Eq. (63) é então imediata 
vindos. 40) 


n. (Z) 


e 
e 


n. (0) 


dd (66) 


Ee o RR À ml 
|+(h/2) 


Um caso particular interessante é h = 2; a atenuação a 
é então inversamente proporcional a n. e o perfil axial 
de densidade é linear. 

A Eq. (66) permite ainda obter uma expressão apro- 
ximada para o comprimento da coluna L, visto que 
n. (L) deve aproximadamente corresponder à densidade 
de corte, n,.. Obtém-se facilmente 


z'iem] 


Fig. 19 Comparação das distribuições axiais calcula- 
das (a cheio) e medidas (a tracejado) da densidade elec- 
trónica numa coluna de plasma de argon mantido por 
O. S. TM, para f = 2450 MHz (1) e 210 MHz (ID e para 
vários valores da potência incidente, em Watts, e da pres- 
são, em Torrs (A - 0,5; B - 0,2; € - 0,1). Em abcissa repre- 
senta-se a distância medida a partir do fim da coluna 
observado experimentalmente 
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l+(h/2) 1 1 + 
- isa 


m (67) 


Ls o 


em que Y,=(w (0)/w)”. 

Apesar das aproximações que foram feitas este mo- 
delo simples descreve com rigor razoável a distribuição 
axial n. (z) e a influência dos vários parâmetros da des- 
carga (N, w, a, P,, etc.) sobre esta distribuição e sobre 
o comprimento da coluna. A Fig. 19 mostra uma com- 
paração dos resultados teóricos com os experimentais 
para diversos valores da frequência, pressão do gás e po- 
tência incidente (neste caso a=0,45cm, b=0,60cm e 
€, = 4.52). Como era de esperar surgem algumas discre- 
pâncias no fim da coluna, especialmente para a frequên- 
cia mais baixa (210 MHz), em virtude das aproximações 
do modelo analítico. Tais discrepâncias poderiam, no en- 
tanto, ser eliminadas se se tivesse em conta o comporta- 
mento correcto da atenuação junto do corte. 


4. CONCLUSÃO 


O presente estado de desenvolvimento da teoria dos 
plasmas produzidos por ondas de superfície permite com- 
preender os processos físicos fundamentais que entram 
em jogo e descrever com bastante rigor as propriedades 
radiais e axiais deste tipo de plasma em função dos pa- 
râmetros da descarga. Neste sentido, o estudo destes plas- 
mas tem permitido progressos significativos na compreen- 
são dos mecanismos essenciais que ocorrem em qualquer 
descarga de radiofrequência ou de micro-onda, abrindo 
desta forma o caminho para a formulação de modelos 
rigorosos e finos de plasmas deste tipo. 

Esta possibilidade revela-se de grande interesse numa 
altura em que as descargas de radiofrequência e de mi- 
cro-onda se começam a usar correntemente no mundo 
tecnológico para a gravura em microelectrónica, a de- 
posição de filmes em superfícies, a excitação de Lasers, 
etc. 
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USO DE MODELOS DE MARKOV 
NÃO OBSERVAVEIS NO RECONHECIMENTO 
DE PALAVRAS ISOLADAS 


António Joaquim Serralheiro”, José Manuel S. Tribolet** 


RESUMO: 


O reconhecimento de fala teve, em especial, a partir 
da última decada um desenvolvimento notável. A finali- 
dade do reconhecimento de fala é a de permitir a um uti- 
lizador dialogar com um sistema não humano. 

Se bem que este objectivo ainda não tenha sido ple- 
namente realizado, os resultados obtidos são encorajado- 
res. O reconhecimento de palavras isoladas pode ser en- 
carado como uma primeira etapa num dispositivo de 
reconhecimento de fala, uma vez que os objectivos a atih- 
gir são mais limitados. 

Este trabalho descreve os sistemas geralmente utiliza- 
dos: o Alinhamento Temporal Dinâmico, DTW, e Mode- 
los de Markov Não Observáveis, HMM, com especial ên- 
fase para este último. 


ABSTRACT 


In the last ten years, speech recognition has greatly 
improved. The final goal for a recognition system is the 
possibility of communication between man and machine. 

Trully, such goal has not been achieved, but results 
are promising. Being less ambicious, isolated word recog- 
nition can be viewed as a first step towards the final 
system. 

This paper describes the two most commonly used 
isolated word recognition systems: DTW (for Dynamic 
Time Warping) and, with special emphasis, HMM (for 
Hidden Markov Models). 


* Assistente do IST. 
** Professor Catedrático do IST. 
Original recebido para publicação em 5/12/86. 
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0. PREÂMBULO 


Este artigo decorre do trabalho de investigação efec- 
tuado pelos autores no âmbito do reconhecimento de 
fala. Nele se fará uma abordagem aos modernos siste- 
mas de reconhecimento acústico global de palavras iso- 
ladas, detalhando-se os princípios teóricos e analíticos 
subjacentes ao uso de modelos de Markov não obser- 
váveis. 

A investigação em Portugal neste domínio é extre- 
mamente recente. Com o objectivo de perspectivar este 
domínio de actividade, um dos autores (J. Tribolet) de- 
dicou um semestre sabático nos Laboratórios Bell nos 
EUA à familiarização dos conceitos envolvidos em re- 
conhecimento de fala e a investigar, em particular, o uso 
de algoritmos de programação dinâmica para alinha- 
mento temporal de padrões (DTW — Dynamic Time 
Warping)!?, cuja tecnologia se encontra hoje amadure- 
cida e já está comercialmente disponível. 

O aparecimento mais recente de sistemas de reconhe- 
cimento baseados em modelos de Markov não observá- 
veis (HMM — Hidden Markov Models) motivou o es- 
tudo e a familiarização de novos conceitos pelos autores, 
os quais são apresentados no presente trabalho. 

Presentemente encontra-se em curso a tese de dou- 
toramento de A. Serralheiro neste dominio, a decorrer 
no INESC, sob a orientação de J. Tribolet e em colabo- 
ração com os Laboratórios Bell. 


1. INTRODUÇÃO 


Neste trabalho far-se-á uma abordagem aos moder- 
nos sistemas de reconhecimento de fala, detalhando para 
o reconhecimento de palavras isoladas. Dos métodos ge- 
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ralmente utilizados, alinhamento temporal de padrões 
(DTW — Dynamic Time Warping) e modelos de Mar- 
kov não observáveis (HMM — Hidden Markov Models), 
referiremos este último mais detalhadamente. Este tra- 
balho principia por situar os sistemas de reconhecimento 
de fala no contexto do processamento digital, aborda al- 
gumas das dificuldades de implementação e seguida- 
mente apresenta o modelo canónico do reconhecedor, 
após o que se pormenorizará cada um dos seus diferen- 
tes blocos operacionais. 

O reconhecimento de fala exige duas fases distintas: 
a aprendizagem ou treino, e o teste ou classificação; na 
primeira o sistema é levado a criar as suas referências 
internas pelo que o(s) utilizador(es) deverá proferir a pa- 
lavra ou palavras que deseje virem a ser reconhecidas; 
a duração de aprendizagem está relacionada com o nú- 
mero de palavras que se quer reconhecer e ainda com 
o desempenho que se pretende para o sistema reconhe- 
cedor. Uma vez concluída esta primeira fase, o sistema 
de reconhecimento encontra-se apto a executar a segunda 
fase que consta do reconhecimento propriamente dito. 
O reconhecimento de fala pode ser encarado ou como 
um problema clássico de reconhecimento de padrões, ou 
através de uma análise semântica do sinal acústico; esta 
segunda versão é por vezes acrescentada à primeira de 
modo a aumentar significativamente a qualidade do sis- 
tema de reconhecimento. 

Este trabalho encontra-se estruturado do seguinte 
modo: no capítulo 2 são apresentadas algumas das mo- 
tivações que impulsionaram o reconhecimento de fala, 
uma breve resenha dos trabalhos pioneiros sobre este as- 
sunto, alguns dos problemas a enfrentar quando se pre- 
tende implementar um sistema de reconhecimento. Refe- 
rem-se os diferentes tipos de reconhecedores, e o modelo 
canónico de um sistema de reconhecimento de fala. 
Ainda sobre este modelo são descritos sumariamente al- 
guns dos blocos constituintes. O capítulo 3 faz referên- 
cia sucinta aos algoritmos de alinhamento temporal di- 
nâmico. O capítulo 4 descreve o recurso aos HMM para 
o reconhecimento de palavras isoladas. O capítulo 5 é 
destinado às conclusões. 


2. BREVES CONSIDERAÇÕES SOBRE 
RECONHECIMENTO DE FALA 


Depois do aparecimento dos computadores digitais, 
começaram a ser desenvolvidos algoritmos de processa- 
mento de sinais, em particular de fala. Com este desen- 
volvimento, tornaram-se possíveis aplicações no campo 
de sintese e do reconhecimento de voz. A possibilidade 
de comunicar oralmente com um computador, ou com 
qualquer dispositivo, começou a ganhar viabilidade, 
sendo hoje comuns sistemas automáticos de resposta e 
de reconhecimento de palavras isoladas. 

Tais sistemas permitem, por exemplo, uma maior ra- 
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pidez de utilização de terminais de computador, uma vez 
que são eliminados os erros de escrita. Deste modo, pos- 
sibilita-se o emprego das mãos para outras tarefas, o que 
é especialmente vantajoso no uso de determinado tipo 
de maquinaria, em que as ordens são dadas oralmente 
sem necessidade de deixar de se prestar atenção ao tra- 
balho em curso. São ainda previsíveis outras aplicações, 
nomeadamente na prestação de serviços, tais como sis- 
temas de reserva de passagens de avião, etc. 

E interessante referir) que o primeiro reconhecedor 
de fala conhecido era um brinquedo — um cão que sal- 
tava quando o seu nome era proferido: “Radio Rex”: 
todavia, quando outros sons vagamente parecidos com 
“Rex” eram proferidos o dispositivo era activado. Pro- 
vavelmente, a primeira tentativa séria de abordar o re- 
conhecimento de fala foi efectuado em 1950 por Dreyfus- 
-Graf em França; o primeiro sistema de reconhecimento 
de palavras isoladas mas que era optimizado para um 
dado utilizador foi construído nos laboratórios Bell em 
1952 por Davis, Biddulph e Balashek permitindo um 
reconhecimento correcto de uma de dez palavras arma- 
zenadas em 97/100 das vezes. O primeiro trabalho que 
utilizava um computador digital data de 1959-1960, rea- 
lizado por Denes e Mathews introduzindo o conceito do 
alinhamento temporal. Hoje, o reconhecimento de fala 
está intrinsecamente ligado ao processamento digital de 
sinais, efectuado quer em computadores de uso genérico 
“main frame general purpose”, quer em processadores 
especializados, nomeadamente microprocessadores, tendo 
aparecido nos últimos anos integrados de elevada capa- 
cidade de processamento que permitem operar em tempo 
real. 

Um sistema de reconhecimento de fala pode ser re- 
presentado pelo seguinte diagrama de blocos: 


») g 


Fig. 1 Sistema de reconhecimento, diagrama de blocos 


wli) 


AQUISIÇÃO 
DE 
SINAL 


PROCESSAMENTO 


As ondas sonoras são captadas por um microfone, 
e o sinal resultante é primeiramente limitado em frequên- 
cia, seguidamente sujeito a um processo de amostragem 
e quantificação, obtendo-se a sequência s(n). Esta sequên- 
cia é, em seguida, processada sendo w(i) o resultado deste 
processamento. 


2.1 Dificuldades 
Sendo a entrada do sistema um sinal acústico, deve- 


rão ser observados certos cuidados na sua recolha, refe- 
rimo-nos à existência de ruído ambiente, por exemplo. 
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Assim, o microfone deverá ser de qualidade, permitindo 
a sua colocação o mais perto possível dos lábios, de 
modo a minimizar a captação do ruído, mas não tão pró- 
ximo para não importunar o utilizador, impedindo a sua 
mobilidade e o seu conforto. No caso de se utilizar o 
micro-auscultador telefónico como entrada de coman- 
dos, é de esperar uma degração do nível de reconheci- 
mento por parte do sistema, uma vez que à baixa quali- 
dade destes transdutores se associam os diferentes tipos 
de ruído presentes nas linhas telefónicas comutadas. 


Mencionámos alguns dos problemas mais óbvios re- 
lacionados com a recolha da informação oral; contudo, 
não podemos deixar de referir que o facto de se preten- 
der alargar o uso destes dispositivos a um maior número 
de individuos, traz problemas adicionais relacionados 
com as variações de dicção existentes entre os diferentes 
utilizadores. Sendo de um modo geral um processo com- 
plexo, o reconhecimento de fala necessita, para ser im- 
plementado, de resposta a um determinado número de 
questõest!): 


|. Que tipo de fala? (continua, palavras isoladas, 
etc.) 

2. Quantos oradores? (um só, alguns, população em 
geral) 

3. Que tipo de oradores? (cooperantes, casuais, trei- 
nados, homens, mulheres, etc.) 

4. Que tipo de ambiente acústico? (sala sossegada, 
ruidosa, lugar público) 

5. Que tipo de equipamento? (microfone de quali- 
dade, micro-auscultador telefónico) 

6. Que necessidade de treino tem o sistema? (algu- 
mas palavras, parágrafos, frases, etc.) 

7. Que treino para os utilizadores? (adaptação na- 
tural, elaborada) 

8. Que dimensão para o vocabulário? (50, 200, 1000, 
10000, palavras seleccionadas, ilimitado) 

9. Que tipo de linguagem? (frases fixas, linguagem 
artificial, português natural, adaptado ao orador) 

10. Que serviço deve ser executado? (respostas fixas, 
etc.) 

1. O que se sabe acerca do orador? (motivação, etc.) 

12. Qual o grau de sofisticação do diálogo” (resposta 
apenas, pedido de repetições) 

13. Que erros se podem tolerar? (nenhuns [<0,1%], 
alguns [<10%], muitos [>20%]) 

14. Qual a rapidez do sistema? (tempo diferido, tempo 
proporcional à frase, tempo real) 

15. Qual o processamento envolvido? (1 Mips, 10 
Mips>10 Mips) 

16. Qual a memória necessária? (1 M, 10 M, + [em 
bit/segundo de fala]) 

17. Qual a sofistificação do software? (programação 
simples, processos concorrentes, processamento 
paralelo, etc.) 


18. Qual o custo por segundo de fala? (1, 10, 100, + 
em milhares de $US) 
19. Data de operacionalidade? 


2.2 Implementação 


Os dispositivos de reconhecimento de fala podem 
agrupar-se, e sem detrimento de outras classificações, 
quanto ao tipo de operação em: 


reconhecedores de palavras isoladas, 
palavras ligadas, 
fala continua, 


e quanto à sua relação com o utilizador em: 


reconhecedores de fala independentes do utilizador, 
dependentes do utilizador. 


Mencionam-se, na literatura da especialidade, dife- 
rentes tipos de reconhecedores, que podemos agrupar em: 


a) reconhecedores de palavras isoladas, em que cada 
palavra é pronunciada isoladamente das outras por 
meio de pausas (silêncios). Temos como exemplo 
“Pára!”, 

b) reconhecedores de palavras ligadas, em que as fron- 
teiras entre as palavras são menos nítidas, como 
por exemplo: *““Três mil e quarenta e um” por opo- 
sição aos anteriores que obrigariam a *'*Três-zero- 
-quatro-um”. 

c) busca de palavras chave, em que se procura de- 
tectar a ocorrência de determinadas palavras, tais 
como: “Ataque”. 

d) reconhecedores de fala contínua, mas num domií- 
nio restrito, em que se tem em vista um determi- 
nado objectivo: “2 lugares para 0 voo...”. 

e) idem, mas qualquer que seja o contexto de apli- 
cação: '*Qual o critério de escolha para adqui- 
Eat 

Í) reconhecimento do utilizador, em que se pretende 
que haja um acesso seguro ao sistema, por ex. apli- 
cavel no “cockpit” de um avião. 


Um sistema de reconhecimento pode ser expresso 
através de um modelo canónico!?) que se indica na 
fig. 2b. 

A entrada do sistema é o sinal acústico proveniente 
de uma fonação (geralmente uma palavra, uma frase), 
fig. 2a. A saída do modelo é a melhor ou as melhores 
estimativas da palavra ou das palavras proferidas. Por 
vezes a saída é ligada a um pós-processador como já se 
referiu. 

O modelo apresentado tem sido largamente utilizado 
devido a várias razões, entre as quais poderemos referir 
que é invariante para vários vocabulários, conjuntos de 
parâmetros, algoritmos de semelhança entre padrões e 
regras de decisão. 
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x(1) 


yCt) yén) 


i Filtro. 


Passa- Baixo 


s(n) 


(8) w 


“À Leisde e 
| decisão 


Extração de | 
parâmetros 


Modelos de | 
| referência 


b) 


Fig. 2 a) aquisição do sinal. b) Modelo canónico para 
reconhecimento de fala 


Por extracção de parâmetros entende-se um processo 
de redução de informação, em que a partir de um con- 
junto extenso de dados, nomeadamente amostras quan- 
tificadas da voz s(n) a um ritmo que poderá ser de 8000 
amostras/segundo (p. ex.), se obtém outro conjunto mais 
pequeno de dados (valores típicos de redução 10 a 100) 
denominados parâmetros, e que descrevem algumas das 
propriedades do sinal acústico. No caso em questão, é 
vulgar o uso de: 


a) espectro localizado (transformada de Fourier), 

b) bancos de filtros, 

c) coeficientes de predição linear, coeficientes Parcor 
(algoritmos de predição linear), 

d) coeficiente cepstrais (modelos homomórficos). 


Optar por qualquer destes, obedece a critérios que 
ponderam: 


a) tempo de cálculo, 

b) memória necessária, 

c) facilidade de implementação, 
d) qualidade a exigir ao sistema. 


Nos sistemas comercializados são vulgarmente em- 
pregues a predição linear e a análise por banco de fil- 
tros, sendo esta última divulgada entre os primeiros sis- 
temas por não necessitar de processamento digital para 
a extracção de parâmetros. 

As regras de decisão permitem, uma vez conhecida 
uma lista de candidatos para a ou as palavras proferi- 
das, obter o melhor candidato segundo critérios pré-es- 
tabelecidost. 3, 


3. ALINHAMENTO TEMPORAL DINÂMICO 


O alinhamento temporal dinâmico (DTW) consiste 
em comparar de algum modo um padrão de teste com 


padrões de referência. Por padrão entende-se uma co- 
lecção de conjuntos de parâmetros extraídos da fala. Por 
sua vez cada conjunto de parâmetros pode representar 
a característica espectral de um segmento de fala. As- 
sim, poderemos dizer que na saída do primeiro bloco 
— extracção de parâmetros, teremos uma sucessão tem- 
poral de “espectros” que de algum modo caracterizam 
a palavra proferida. Como referências teremos, e para 
cada palavra do vocabulário, padrões que foram obti- 
dos através de um processo de aprendizagem. 

Consiste, este método, em determinar a “semelhança” 
entre os padrões de teste e de referência, o que é efec- 
tuado através de uma distância. Como a duração dos 
padrões varia, torna-se necessário um alinhamento tem- 
poral destes antes de se proceder à determinação da se- 
melhança entre os mesmos. Consideremos dois padrões, 
um de teste T(t) e outro de referência R (1); o objectivo 
é determinar a função de alinhamento W (t), ou de “'war- 
ping”, que efectua a correspondência de R para T, ba- 
seando-se na minimização da distância entre os dois pa- 
drões: 


D (T,R) = min / dit, W(D] G [t, W(D, W'(D] dt 
tw (0) 


sendo [W (t)l o conjunto de todas as funções monoto- 
nicamente crescentes, com derivada contínua, d[t, W (t] 
uma distância pontual de R a T, G é uma função de pon- 
deração. 

Seja T=/(T(1), T(2),..., T(NT)) o conjunto dos pa- 
drões de teste e R= [R(1), R(2),..., R(NR)! o conjunto 
dos padrões de referência. O alinhamento temporal óp- 
timo é a curva que relaciona o eixo temporal m do pa- 
drão de teste T (m) com o eixo dos tempos n do padrão 
de referência, na forma: 


m=W(n) 


Como cada padrão é constituído por um conjunto 
contável de “espectros”: tantos quantos os segmentos 
em que se partiu o sinal correspondente a cada palavra, 
a determinação da função de alinhamento consiste na 
determinação de um caminho numa rede. Esta determi- 
nação envolve a especificação de vários factores: 


a) restrições aos pontos inicial e terminal, 
b) restrições locais: “percursos possíveis”, 
c) restrições globais: limitação do espaço, 
d) escolha dos eixos para os padrões, 

e) medida das distâncias. 


Não detalharemos os pontos anteriores por se encon- 
trar fora do âmbito deste trabalho, veja-se! 3 19, 
3.1] Implementação do algoritmo de DTW 


Para se implementar o algoritmo de alinhamento tem- 
poral dinâmico, utilizando um procedimento de progra- 
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mação dinâmica, consideram-se válidas as seguintes as- 
serções: 
a) um percurso globalmente óptimo também o é lo- 
calmente, 
b) o caminho óptimo para o ponto da rede (n,m) só 
depende dos valores de n', m' tais que: nº <n, 
m'<m. 


Define-se assim a relação recursiva de programação 
dinâmica; seja D, (nm, m) a distância parcial acumulada 
desde (1,1) até (n, m): 


D, (n, m=d(T, R9)+min [D, (n-1, q) 
q<m 


D, (n, m) pode ser definida recursivamente em ter- 
mos de um ponto intermédio n', m' [n'<n, m'<m]: 


D, (n,m)=min [D,(n', m”)+ d((m', m),(n, m))] 
(n';m') 


em que d é a distância ponderada de (n”, m')-(n,m) e 
L é o número de segmentos no percurso anterior, e W (n) 
é uma função de ponderação: 


à Um), (nm) =D di(n-D,(m-DIW (nº) 


Para uma implementação eficiente de D, (n,m) é 
apenas necessário aplicar restrições à variação do per- 
curso de (n',m”) para (n,m) dentro da rede. Uma das 
restrições utilizadas é a de Itakura, em que os percursos 
possíveis de (n',m') para (n, m) são os indicados abaixo; 
note-se que o ramo assinalado com x corresponde a não 
permitir que o percurso seja horizontal durante mais do 
que um segmento. 


Fig. 3 Percursos possíveis, restrição de ITAKURA. 


D, (mm) = (min/D, (n-1, m-2)+ d(n,m), 
D,(n-1, m-D+d(n,m),D,(n-1,m)g(k)+ d(n, m)) 


sendo g (k) = «o no ramo assinalado com x, caso con- 
trário, g (k) = 1. 

Para a implementação do algoritmo teremos: 

a) início: DA, D=d(,1) W(1) 

b) recursividade: calcule D, (n,m) recursivamente 

para I<n<N, I<m<M 
c) finalização: 
l 


D=D.(NM=—— 
a É N(W) 


da 


sendo N(W) uma função de normalização. 


3.2 Algumas considerações sobre o DTW 


Para ilustrar a determinação de D, (n,m), atente-se 
no seguinte exemplo, em que temos para a palavra de 
teste FOUR pronunciada FORR e que como referência 
temos FFOOR, e admitamos que cada letra é um bloco, 
que os pontos terminais foram correctamente determi- 
nados, e que se obtinha a seguinte matriz de verosseme- 
lhança: 


R|9 622 
. 0/8186 
Zols18s 
SF|2377 
e ÉiiTAs 

FORR 
(Testa) 


Fig. 4a) Matriz de verossemelhança. 


da qual se obtém a matriz das distâncias acumuladas mí- 
nimas: 


Rizs 11.7... 
= 0 |19 5 13 18 
& 0/11 4 12 18 
S F|3 9 15 22 
O ' 
“ F|1 815 20 


IF O RR 


(Teste) 


Fig. 4b) Matriz das distâncias acumuladas mínimas 


Teriamos deste modo que a distância acumulada mi- 
nima entre as duas palavras é de 9. Procederíamos de 
igual modo para os restantes padrões de referência, ob- 
tendo várias distâncias. Estas seriam sujeitas a um pro- 
cesso de decisão que elegeria um ou mais candidatos para 
a palavra FORR, e se não ocorresse nenhum erro de re- 
conhecimento, o sistema reconhecê-la-ia como FOUR. 

O desempenho de sistemas de reconhecimento de pa- 
lavras isoladas é bastante bom, havendo sistemas comer- 
cializados em que a taxa de reconhecimento correcto 
atinge os 99%, para um só utilizador. O aparecimento 
de circuitos integrados dedicados ao reconhecimento por 
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DTWºC! vem permitir uma maior difusão dos sistemas 
de reconhecimento. 

Quando a recolha do sinal acústico é efectuada num 
ambiente ruidoso, quando se trata do reconhecimento 
de palavras ligadas, ou ainda com alguns sons pode 
tornar-se difícil a localização exacta dos início e fim das 
palavras, e como tal podem ocorrer erros no reconheci- 
mento, uma vez que se estão a atribuir distâncias eleva- 
das aos primeiros segmentos dos padrões, e eventual- 
mente podem ser logo excluídas as palavras correctas. 
veja-se, p. ex.) para uma abordagem mais completa. 


4. OS MODELOS DE MARKOV NÃO-OBSERVÁVEIS 
E O RECONHECIMENTO DE PALAVRAS 
ISOLADAS 


Os processos de reconhecimento acústico global têm 
como base modelos de produção de fala; através de ob- 
servações efectuadas na saída do aparelho fonador tenta- 
-se inferir qual a palavra ou conjunto de palavras pro- 
feridas, pertencentes a um dado vocabulário e que são 
o suporte da informação que se pretende transmitir. 

De uma forma necessariamente simplista podemos 
supor que as palavras são proferidas a partir de um nú- 
mero finito de configurações do tracto vocal, e é o modo 
como essas configurações se sucedem ao longo do tempo 
que produzem os diferentes sons. Assim sendo, podere- 
mos ainda supor que se num dado instante o tracto vo- 
cal está numa dada configuração há uma certa proba- 
bilidade de no instante de tempo seguinte se encontrar 
noutra configuração. 

A ocorrência dos sons de uma dada palavra está, por- 
tanto, condicionada a um processo probabilístico. So- 
mos deste modo levados a admitir um modelo probabi- 
lístico para o processo de produção da fala. A variação 
temporal das caracteristicas espectrais do sinal de fala 
pode ser associada à lei probabilística de transição de 
estados de uma cadeia de Markov, enquanto que para 
cada instante de tempo, o espectro observado é condi- 
cionado a uma observação efectuada num dado estado 
ocupado. A ordem do modelo de Markov poderá ser 
qualquer, e as saídas dos seus estados podem ser um pro- 
cesso aleatório multivariavel com densidades de proba- 
bilidade conjunta contínuas. 

No que se segue restringir-nos-emos a modelos de pri- 
meira ordem e a observações que foram retiradas de um 
alfabeto discreto e finito de acordo com funções de dis- 
tribuição discretas associadas a cada estado. 

Começar-se-á por uma referência a conceitos funda- 
mentais: 

Definição: Um processo de Markov de 12 ordem é 
um processo estocástico cujo passado não tem influên- 
cia no futuro se o presente é especificado”: 


se t,.j<t, então 
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PX()<x | X(O, t< tr) =PXO)<x | X(t.) 


donde se retira que se t;<t,<...<t, 
teremos: 


PX(O)<x | XL Does XD) =PXC()<x, | X(t,.,)) 
e para sistemas discretos no tempo: 
PIX, <xa | Xacares Xi) = PÍX, x, | Xe] 


Definição: Uma função probabilística de uma cadeia 
de Markov não observável é um processo estocástico ge- 
rado por dois mecanismos interligados!!?: 


— uma cadeia de Markov subjacente com um nú- 
mero contável de estados; 

— um conjunto de funções aleatórias, cada uma das 
quais está associada a cada estado. 


Assim, num dado instante n o processo é suposto es- 
tar num só estado e um símbolo é produzido de acordo 
com a função de distribuição de probabilidade que de- 
pende do estado actual. A cadeia de Markov subjacente 
muda de estado de acordo com a matriz probabilística 
de transição. O observador apenas vê como saídas (os 
simbolos anteriormente referidos) valores das funções 
aleatórias associadas a cada estado e não pode obser- 
var directamente os estados da cadeia de Markov, donde 
o nome de ''não observável”. Os valores das funções 
aleatórias correspondem, por exemplo, a conjuntos de 
coeficientes de LPC, associados a uma dada configura- 
ção do tracto vocal. 

Suponha-se que uma cadeia de Markov tem N esta- 
dos q,, q, q, € que as observações são retiradas de 
um alfabeto V de M espectros v,, V,,., Yy: À cadeia 
pode ser especificada em termos de um vector de distri- 
buição de estados inicial w' =(7,, Ts» Tx) € UMa ma- 
triz de transição de estados A =[a,] I<i, j<N. O sig- 
nificado de 7, é o da probabilidade de nos encontrar- 
mos num estado q, num instante de tempo inicial, t = 0. 
Quanto a a, j é a probabilidade de transitar (no instante 
seguinte) para o estado q, estando no estado q;: 


a;=P(qg, em t+1 | q, em t) 


QREN 


Fig. 5 Probabilidades de transição 


O processo aleatório associado com os estados pode 
ser representado por uma matriz estocástica B= Ib,d 
com ISj<N I<Sk<M sendo b,, a probabilidade de 
observar o símbolo v, sabendo que é q, o estado actual: 


by =P(v, em t|q, em t) 
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Deste modo, o HMM é caracterizado pelo terno or- 
denado (7m, A, B). 

Para se poder utilizar os HMM para reconhecimento 
é necessário resolver dois problemas: 


— estimação da sequência de observações, que será 
usada para a classificação de uma dada fonação, 

— estimação dos parâmetros do modelo e que ser- 
virá para o treino deste. 


Ambos os problemas provêm de uma sequência de 
observações O, O,... O, em que cada O, (I<t<T) é al- 
gum v, e V, isto é uma observação correspondente a 
um símbolo do vocabulário. 

Quanto à classificação, pretende-se reconhecer pala- 
vras de um vocabulário seleccionado de W palavras 
W,, Wo;,.., W,. Como entrada do sistema tem-se uma 
sequência de observações O e um conjunto de v mode- 
los (um por cada palavra, pelo menos) M,, M,,..., M,. 
Tem de se calcular a probabilidade de a sequência ser 
gerada por cada um dos modelos: P,=P(0/M,) 
(I<i<V). SeP,>P, I<xj<Va palavra de teste será re- 
conhecida como a palavra W,; o mesmo é dizer que se 
reconhece como palavra pertencente ao vocabulário 
aquela que tiver a maior probabilidade de ter sido ge- 
rada pelo modelo i. 

Quanto ao treino do sistema de reconhecimento, 
trata-se de determinar os modelos M,=(r,, A, B, 
(I<i<V), a partir de K referências de treino: OM, 
092,.., O*, sabendo que a observação O? é prove- 
niente de uma palavra (falada) W, do vocabulário. 


4.1 Classificação (teste) 


Uma vez que o HMM lida com um número finito 
de simbolos, para se poder reconhecer fala recorrendo 
a este método torna-se necessário obter um conjunto nu- 
merável de símbolos. Assim, e se se estiver a empregar 
coeficientes de predição linear, estes terão de ser sujei- 
tos a um processo de quantificação, por exemplo, a re- 
ferida quantificação vectorial). (1), Recorre-se à quan- 
tificação vectorial em vez de quantificação escalar, 
porque para uma mesma distorção aquela necessita de 
um menor número de índices. 

Seja O=0, O,... O, uma sequência finita de obser- 
vações onde cada observação é um simbolo discreto re- 
tirado de um alfabeto finito de símbolos: os índices de 
quantificação vectorial. 

Dada uma sequência de observações O e o modelo 
M = (7, A, B) com N estados, a probabilidade de O ter 
sido gerada por este modelo é!M; 

P(O | M) “+ Ty Di (0) Aee Ary Dip (Or) 


01 EF 


O significado da equação anterior é o seguinte: Ini- 
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cialmente, t= o, estar-se-á no estado i, com probabili- 
dade 7,,, e é “gerado” o símbolo o, com probabilidade 
b, (0). Em seguida transitar-se-á para o estado se- 
guinte, i,, com probabilidade a, ,., e gera-se o simbolo 
o, com probabilidade b, (o,). Este processo prosseguirá 
até se atingir o último estado. 

O somatório anterior pode ser calculado de um modo 
mais eficiente se definirmos probabilidade progressiva, 
a, (i): 


a,(i) = P(o, 0,... O, q em t|M) 


isto é, a probabilidade de observar a sequência o, 0,... 
0,.,, estando no estado q, num dado instante t, sa- 
bendo ou dado um modelo M. 

Pode-se, agora, obter a seguinte relação de recursi- 
vidade: 


N 
D+ l ()) = IL a, (i) a,] b, (o,, ) t == l, 2, nu. 


Definindo probabilidade regressiva: 


E D=(6,,,0, or |q, em t, M) 
como sendo a probabilidade de observar a sequência 
0,,1 9,,2:- Oy Sabendo que se ocupava o estado q, em 
t e um dado modelo M, obteremos a seguinte relação 
de recursividade: 


N 
8 D=L ab; (0,1) Ba O) T-1>t>1 


e poderemos então determinar a probabilidade de P na 
forma: 
mM N 
P=P(O |] M)=L, dr q, (1) ab; (0,,,)) Br, O) 
I<t<T-1 


ou ainda, utilizando apenas probabilidades progressivas: 
N | 
P= L a (i) 


Ao utilizar a equação anterior para o cálculo de P 
está-se a incluir implicitamente todas as sequências de 
estados dentro de um determinado modelo que pode- 
riam ter gerado a observação. 


Poder-se-ia definir P como o máximo, para todas as 
sequências de estados i=iç, à,» ày da probabilidade 
conjunta P(o, i). Esta sequência de estados e a corres- 
pondente probabilidade da sequência de observação po- 
dem ser determinadas recorrendo ao algoritmo de Vi- 
terbitD:; 

início: ISISN 


6, ()=7,b, (0,) 
V, (D)D=0 
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recursão: à, (j) = max [ó.., (i) a] b, (0) 
V, ()=arg max [ó,., (1) aj] 


2<t<T ISj<N 

Este algoritmo, além de permitir determinar qual o 
modelo que poderia ter dado O como observação, de- 
termina ainda qual a melhor sequência de estados que 
origina P,.. 


finalização: P* = max [ó, (i)] 
ISI<N 


E a : 
Hg = arg Toa [dr ()] 


Para se obter a sequência de estados mais provável: 
para t=T-1, T-2,.., 1: 


= Wi Wy) 


Repare-se que se está a caminhar do último estado 
para o primeiro estado da cadeia (''backtracking””). 


4.2 Treino do modelo 


Sendo impossível obter um conjunto infinito de ob- 
servações, o modelo que se irá obter não será óptimo 
pelo que normalmente, o modelo obtido corresponderá 
a um máximo local de P. Consideremos as probabilida- 
des progressiva e regressiva, e determine-se o número ex- 
pectável de transições y, de q, para q,, condicionada à 
sequência observada: 


<A 
hp 


T-1 
o q (1) a db; (0,9) Bi, QU) 

Deste modo, o número expectável de transições do 
estado q, é dado por y;: 


«Lt mês 
n=p E, % (1) 6, (1) 


— 


N 
Recorde-se que %, == Yi 


O quociente 7,;/y, representa uma estimativa da pro- 
babilidade de nos encontrarmos no estado q, sabendo 
que o anterior era q,. Este quociente é uma estimativa 
da probabilidade de transição a, 


T. 
= Ay (1) dj b, (0, )B,.1 0) 

O CR o 
Yi bo C, (1) B, (i) 


De modo semelhante se pode determinar uma esti- 
mativa de b,, como a frequência de ocorrência de v, no 
estado q;; teremos, em termos de aq e de 5: 


L a a, (4) 6, Q) 
EE Es 


— 130,= 


DRM ONO 


ik 


E, finalmente, para a probabilidade inicial II: 


= « (i) 8, (1) 


l 
P 

Demonstra-se!'?) que estas estimativas incrementam 
P excepto num ponto crítico que corresponde a máxi- 


mos locais. As equações anteriores derivam do algoritmo 
de Baum-Welch (B.W.)t2, 


4.3 Estruturas de HMM e número de estados 
Consideram-se, pela sua generalidade, três tipos de 


estruturas: não restringidas e restringidas em cascata e 
paralelo, fig. 6: 


Fig. 6a) Estruturas de cadeias de Markov, 


não restringidas 


Fig. 6b) Estruturas de cadeias de Markov, restringidas, 
em cascata 


ABRIL 87 47 


Estruturas de cadeias de Markov, restringidas, 
em paralelo 


Fig. 60) 


Os modelos com as restrições apresentadas têm ma- 
trizes de transição triangulares superiores; também de- 
nominados modelos esquerda-direita. Os modelos de es- 
trutura em cascata podem progredir sequencialmente da 
esquerda para a direita, sendo possível o salto de um ou 
mais estados; os de estrutura em paralelo permitem ca- 
minhos múltiplos ao longo do modelo, continuando a 
evoluir da esquerda para a direita. 


Qualquer das estruturas anteriores pode apresentar 
um número de estados arbitrário; notemos, contudo, que 
o número de parâmetros evolui com Nº (matriz A), N*M 
(matriz B), e se é utilizado um número de estados grande 
(=10) o cálculo de P vem fortemente penalizado. Repare- 
-se que o modo de teste é o mais utilizado, pelo que os 
cálculos deverão ser rápidos; quanto ao treino do mo- 
delo, pode correr fora do período de utilização do sis- 
tema, e ocorrerá, geralmente, antes de o utilizar pela pri- 
meira vez, e eventualmente, para posteriores actualizações 
do vocabulário. Aparentemente não há um processo sim- 
ples de determinar quantos estados deverá ter o modelo, 
uma vez que os estados de cadeia de Markov não têm 
que estar relacionados fisicamente com o fenómeno em 
observação. 


Rabiner et al!) para o reconhecimento dos dez di- 
gitos determinaram a evolução do erro de detecção para 
um número de estados variável de 2 a 20 em estruturas 
em cascata, fig. 7: 


Fm 


Da 


Erro 


Fig. 7 Evolução do erro com N 
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Nota-se uma diminuição do erro com o aumento do 
número de estados, contudo, esse decremento é menor 
a partir de 8-9 estados, pelo que o acréscimos de com- 
plexidade do modelo não é justificado pelos resultados 
obtidos. As flutuações do erro em torno do valor mé- 
dio são devidas à escolha inicial das matrizes A e B. 

Na figura seguinte, são representadas as evoluções do 
erro de reconhecimento para cada um dos dígitos em fun- 
ção do número de estados, e mais uma vez se verifica 
que para este exemplo concreto, 5 ou 6 é o número de 
estados que optimiza o erro de reconhecimento bem 
como a velocidade de processamento!!», Note-se ainda 
que para palavras monossilábicas (“'one, two””) o erro 
é maior do que para palavras bissilábicas (“'zero, seven”), 
o que leva a crer que não existe relação simples com o 
número de estados e a duração de sons elementares!" 


To pç E E tf=E a 


Fig. 8 Evolução do erro com N para cada um 
dos dígitos 


4.4 Sequências múltiplas de observação 


Numa tentativa de aumentar o desempenho do sis- 
tema de reconhecimento quer para sistemas mono-utili- 
zador que para multi-utilizador, é acessório ter várias ob- 
servações de uma mesma palavra de modo a caracterizar 
as suas diferentes realizações. Tal facto irá penalizar bas- 
tante o processo de aprendizagem. Realizações práticas 


TÉCNICA 3/4-86 


utilizam de 2 a 10 modelos para a mesma palavra a ser 
reconhecida, o que pode implicar um aumento do treino 
de cerca de uma ordem de grandeza. E claro que o teste 
ou classificação também é afectado, contudo, o desem- 
penho superior pode compensar o acréscimo de proces- 
samento, o que dependerá da aplicação em vista. 


4.5 Restrições durante o treino 


Para estruturas não restringidas, os elementos das ma- 
trizes A e B podem assumir quaisquer valores; nas res- 
tringidas, alguns desses elementos serão forçosamente nu- 
los. Assim e para os modelos da fig. 9 teremos: 


Fig. 9 Modelos de Markov restringidos 


Quanto à matriz B, e pelo facto de as estimativas ini- 
ciais poderem originar algum b, (k) = 0, pode suceder 
que «,., (i).a, seja diferente de zero para apenas um va- 
lor de j, e como tal a probabilidade dessa observação 
(0, = k) provir do modelo é nula — o que corresponde 
a um erro de falta de aprendizagem ou aprendizagem 
insuficiente. Um modo de resolver esta questão consiste 
em não anular valores de b, (k) cuja estimativa dê nula, 
ou seja impor b, (k) > e sendo «e um limiar escolhido. 
Teremos que os b, (k) inferiores ao limiar terão o seu 
valor modificado para e, os b, (k) superiores a e são es- 
calados de modo a manter - b, (k) = 1. 


Na figura seguinte, fig. 10 o erro de reconhecimento 
é obtido para dois modelos com estruturas distintas!): 
um é não restringido, e o outro é-o, ambos têm 5 esta- 
dos. A variável independente é o limiar «e que varia de 
O a 10º, verificando-se um erro de reconhecimento 
mais elevado, o que confirma o atrás exposto. Uma vez 
que o modelo não restringido tem um desempenho in- 
ferior serão utilizados apenas modelos restringidos, e des- 
tes os em cascata. Para b, (k) qualquer valor não perto 
de zero é aceitável, e como exemplo, poderemos ter 
b, (k) = 10º. 
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Fig. 10 Erro de reconhecimento, modelos restringido 


e não restringido, 5 estados 


4.6 Escalamento 


O método de aprendizagem apresentado necessita do 
cálculo de a e de 8 para I<t<T, I<i<N, contudo, 
das equações que nos permitem determinar estas pro- 
babilidades verifica-se que elas tendem exponencialmente 
para zero quando T co. Deste modo, para se treinar 
convenientemente um modelo, ocorrerá uma situação de 
“underflow” que pode ser resolvida mediante um esca- 
lamento conveniente de « e de 8. Baseia-se este escala- 
mento no produto de « por um coeficiente independente 
do estado de modo a que esta permaneça dentro da gama 
dinâmica do computador; 8 sofrerá um escalamento 
idêntico. Para uma descrição mais pormenorizada do es- 
calamento das probabilidades progressiva e regressiva, 
veja-set!?, 


4.7 Modelos de Markov esquerda-direita 


O uso de modelos que evoluam sequencialmente da 
esquerda para a direita tem-se revelado muito útil nos 
sistemas de reconhecimento de fala. Apresentam as se- 
guintes propriedades; 


1) a primeira observação é produzida quando o mo- 
delo se encontra num determinado estado, deno- 
minado estado inicial, q,. 

2) a última observação é produzida com o modelo 
num estado final, de onde não transitará a não ser 
para ele próprio (absorvente), q,. 

3) uma vez que um estado seja abandonado, não será 
mais ocupado. 


Nestas condições teremos: 
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Fig. 11 


Modelo de Markov esquerda-direita 


No caso de se utilizar um modelo nestes pressupos- 
tos, para a sua aprendizagem não interessa uma sequência 
de observação muito longa, porque uma vez atingido o 
estado final o que resta da sequência de observações não 
traz nenhuma informação sobre os estados anteriores. 

O conjunto de treino para este tipo de modelos con- 
sistirá num conjunto de observações obtidas por dife- 
rentes começos no estado inicial. Se se tratar do reco- 
nhecimento de palavras isoladas, várias fonações da 
mesma palavra constituirão o conjunto de treino. Deste 
modo, e para calcular a probabilidade de uma sequên- 
cia de observações 0,, 0,,..., Oy ter sido gerada pelo mo- 
delo, sabendo ainda que (e de acordo com as proprie- 
dades 1 a 3) que o, e o, teriam sido produzidas nos 
estados inicial e final, respectivamente: 


[= (1, 0,..., 0) 


B, (D=1 l <j <N 

Quanto à re-estimativa de a, para satisfazer a con- 
dição 3, basta impor a, =0 para j<i e para outras situa- 
ções que se entenda manter (por exemplo, o número má- 
ximo de estados a saltar). Esta condição é suficiente para 
manter as re-estimativas das probabilidades de transição 
nulas se se utilizar o algoritmo de Baum-Welch. Con- 
tudo, este algoritmo terá de ser adaptado às multiplas 
observações da mesma palavra, pelo que teremos: 


Seja O = |O, 082,..., OM] o conjunto de sequências 
de treino em que O%* é a sequência indice k. 


Todas as observações são supostas serem independen- 
tes entre si, e ajustar-se-ão os parâmetros do modelo M 
para maximizar: 


a K 
P=II P(OV|/M)=IH P, 
k=1 k=1 


Uma vez que o algoritmo de BW. determina a fre- 
quência de ocorrência de vários acontecimentos, é ape- 
nas necessário calcular cada uma dessas frequências re- 
lativas em cada sequência separadamente e adicioná-las 
para todas as sequências. Teremos assim: 


k i 
b o 


mm 


2 (1) a; b; (OP) BM, Q) 
K ts | 
E Lo a (i) BL (i) 
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| “a () BM () 
asim 
E Lo at () BM (1) 


Note-se que o vector II, não é reestimado. 


6. CONCLUSÕES 


Como se referiu anteriormente, o algoritmo de BW. 
não garante a obtenção de máximos globais de P(0/M) 
mas apenas de máximos locais. Contudo, e por inspec- 
ção de P(O/M) pode-se determinar, para um modelo 
esquerda-direita, a sequência de estados mais provável 
de comprimento T e que será: 


Sa Sigo (T-N) O Syqueyr+ees Sin) COM T>N 


sendo R O S, a sequência de R repetições do estado i. 
Esta sequência pode ser fornecida como condição ini- 
cial para o treino do modelo e conduzirá com maior pro- 
babilidade a um modelo que estará mais perto do óp- 
timo do que um obtido com uma condição inicial 
aleatória. A escolha de um modelo esquerdo-direita para 
modelo a utilizar tem a sua razão assente quer num me- 
nor número de cálculos necessários para o treino do mo- 
delo (restrições na matriz A), quer numa tentativa de fa- 
zer corresponder a uma evolução temporal da fala numa 
evolução do modelo da esquerda para a direita. 

Os modelos de Markov aplicados ao reconhecimento 
de palavras isoladas têm dado bons resultados quer em 
valor absoluto: desempenhos superiores a 90%, quer 
quando comparados com outros métodos de reconheci- 
mento, nomeadamente o DTW, especialmente se este 
também incluir quantificações vectorial dos parâmetros. 
Se os HMM utilizarem como saidas dos seus estados pro- 
cessos aleatórios multivariáveis com densidades de pro- 
babilidade conjunta continuas, o seu desempenho é idên- 
tico ao dos DTW. A vantagem da utilização dos HMM 
reside ainda numa maior rapidez do processo de classi- 
ficação em que é francamente superior ao DTW; con- 
tudo, o treino é, do ponto de vista do tempo necessário 
mais demorado relativamente ao DTW. Para aplicações 
correntes este limitação não determinante. 

Resta referir alguns dados sobre as prestações típi- 
cas de reconhecimento baseadas no DTW e no HMM: 
assim, e de acordo com resultados obtidos nos Labora- 
tórios Bell, teremos que, para reconhecer um entre os 
dez digitos se obtém uma taxa de sucessos de 98,2% para 
o DTW e de 97,1% para o HMM este último com quan- 
tificação vectorial (dimensão do codebook: 128 code- 
words). Estes valores dizem respeito a sistemas multi-uti- 
lizador. 
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A ESTABILIDADE 
ESTRUTURAL E O 


COMPORTAMENTO 
POS-ENCURVADURA 


A.J. REIS 


Por lapso verificado na montagem do último número da Re- 
vista Técnica, não foi inserido no artigo A ESTABILIDADE ES- 
TRUTURAL E O COMPORTAMENTO PÓS-ENCURVADURA 
da autoria do Prof A. J Reis, o quadro que agora reproduzimos. 

Pelo facto apresentamos as nossas desculpas. 


QUADRO I 


Ponte sobre o Danúbio em Viena 
colapso em 6-NOV.-1969 


Banzo comprimido do caixão — Placa de 7600 x 10mm 
com 12 reforços de 160 x I2mm; 
o. = 244 N/mm? 


TENSÃO MÉDIA DE 


IMPERFEIÇÕES | 
COLAPSO LOCAL 


do (mm) + 0, (N mm?) 


WALKER & 
MURRAY 


HORNE & 
NARAYANAN 


FRIEZE 
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COLECÇÃO DE PROGRAMAS E FILTROS 


DESTINADOS A OPERAR 


SOBRE “CONJUNTOS VAGOS” 


A. J. N. Martins Aguiar* 


RESUMO 


Colecção compreensiva de programas aplicáveis a 
conjuntos vagos (Fuzzy sets) e respectivas demons- 
trações. 


ABSTRACT 


Comprehensive collection of programs applicable 
io Fuzzy sets and respective proofs. 


Apresentação 


O objectivo da apresentação deste programa e 
respectiva justificação, é pôr à disposição do utente 
um conjunto de «filtros» bem como das conectivas 
usadas na Teoria dos Conjuntos Vagos (Fuzzy Sets, 
Ensembles Flous, Conjuntos Borrosos). 

Na verdade, sem estas subrotinas torna-se pesado 
operar no domínio dos conjuntos vagos. 

Este trabalho além de resumir alguns trabalhos 
anteriores realizados nas cadeiras de Automação | 
e II, cujos nomes são referidos na Bibliografia, 
introduz todo um conjunto de demonstrações for- 
mais e de novas subrotinas. 

Toda esta matéria foi programada em Fortran IV 
porque é uma linguagem relativamente geral e dai 
que susceptível de correr em computadores de diver- 
sos tipos, e é assim possível usar este conjunto de 
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programas na resolução de problemas propostos em 
«Fuzzy Mathematics». 

Outros programas destinados a completar as apli- 
cações dos Conjuntos Vagos serão apresentados em 
futuros artigos. 

Os fundamentos teóricos estão referidos na Biblio- 
grafia. 


A. G. PORTELA ** 


Introdução 


As conectivas definidas nos conjuntos vagos são 
apoiadas nas conectivas definidas nos reticulados onde 
se descrevem os respectivos «caracterizantes». 

Assim é essencial por um lado identificar os reti- 
culados referidos para escolher as conectivas aplicá- 
veis, a fim de operar sobre os conjuntos vagos sub- 
jacentes. 

Nesta ordem de ideias será apresentado: 


f.º Capitulo 


CARACTERIZAÇÃO DO CONJUNTO ONDE 
UMA RELAÇÃO DE ORDEM FOI IM. 
POSTA 


1.1] — Introdução. 
1.2 — Verificação das propriedades do grafo 


a) Existência e unicidade de ínfimo e su- 
premo do grafo; 

b) Existência de circulações; 

c) Existência de ramificações; 

d) Resumo e conclusão de 1.2. 
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